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The quest for experimental platforms that allow for the exploration, and even control, of the interplay of low 
dimensionality and frustration is a fundamental challenge in several fields of quantum many-body physics, such 
as quantum magnetism. Here, we propose the use of cold crystals of trapped ions to study a variety of frustrated 
quantum spin ladders. By optimizing the trap geometry, we show how to tailor the low dimensionality of the 
models by changing the number of legs of the ladders. Combined with a method for selectively hiding of ions 
provided by laser addressing, it becomes possible to synthesize stripes of both triangular and Kagome lattices. 
Besides, the degree of frustration of the phonon-mediated spin interactions can be controlled by shaping the trap 
frequencies. We support our theoretical considerations by initial experiments with planar ion crystals, where a 
high and tunable anisotropy of the radial trap frequencies is demonstrated. We take into account an extensive 
list of possible error sources under typical experimental conditions, and describe explicit regimes that guarantee 
the validity of our scheme. 
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I. INTRODUCTION 

The transition from single-particle to many -body quantum 
systems yields one of the most interesting concepts of physics: 
emergence. As emphasized by R W. Anderson [1], the laws 
that describe the collective behavior of interacting many-body 
systems can be fundamentally different from those that gov- 
ern each of the individual particles. This effect leads to some 
of the most exotic phenomena of condensed-matter physics 
in the last decades [2]. Unfortunately, the high complexity 
of many -body systems turns our endeavor to understand such 
emergent phenomena into a fundamental challenge for both 
experimental and theoretical physics. Since exact analytical 
solutions seldom exist, even for oversimplified models, one is 
urged to develop efficient numerical methods. An alternative 
to this approach are the so-called quantum simulations [3], 
which make use of a particular quantum system that can be 
experimentally controlled to a high extent, in order to unveil 
the properties of a complicated interacting many-body model. 

There are two different strategies for the quantum simula- 
tion of many -body systems [4], the so-called analog and digi- 
tal quantum simulators (QSs). On the one hand, an ideal ana- 
log QS would be a dedicated experimental device where one 
can prepare the quantum state, engineer a Hamiltonian of in- 
terest, and measure its properties. On the other hand, a digital 
QS aims at reproducing the dynamics of any given Hamilto- 
nian by concatenating a set of available quantum gates. In- 
dependently of their particular experimental implementation, 
these types of QSs would be capable of exploring models 
from very different areas of physics, ranging from condensed- 
matter to high-energy physics. Let us remark that this enter- 
prise benefits directly from the development of architectures 
for quantum-information processing [5]. In this work, we fo- 
cus on one of these architectures: laser-cooled atomic ions 
confined in radio-frequency traps [6]. Here, some remarkable 
quantum simulations have already been accomplished in ex- 
periments, either in the digital [7] or analog [8-12] versions. 
In particular, we shall concentrate on analog QSs, where a 
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number of theoretical proposals already exist [13]. Some of 
these proposals target the quantum simulation of relativistic 
effects [14], quantum spin models [15-18], many-body boson 
systems [19], spin-boson models [20], or theories of quantum 
transport and friction [21]. 

A direction of research that is being actively explored is the 
QS of magnetism [15], which yields a playground for a vari- 
ety of cooperative quantum many-body effects. Even though 
the phenomenon of magnetic ordering was already known to 
the ancient greeks, the understanding of its microscopic ori- 
gin could only be achieved after the development of quantum 
mechanics. A particular topic of recurring interest is the fate 
of the magnetically-ordered phases in the presence of quan- 
tum fluctuations [22], which get more pronounced as the di- 
mensionality of the system is reduced. Such quantum fluc- 
tuations may destroy the long-range order, favoring paramag- 
netic phases and triggering the so-called quantum phase tran- 
sitions [23]. Alternatively, they may stabilize more interest- 
ing phases, such as the long-sought two-dimensional quan- 
tum spin liquids, which have connections with the theory of 
high-temperature superconductivity [24] . The advent of a QS 
for quantum magnetism would allow for an unprecedented ex- 
perimental realization of these phases, addressing puzzles that 
remain unsolved due to their great complexity. 

A representative model that exemplifies the usefulness of 
analog QS is the one-dimensional quantum Ising model [25], 
which describes a chain of interacting spins subjected to 
a transverse magnetic field. This model has traditionally 
been considered as a cornerstone in the theory of quan- 
tum phase transitions [23]. However, the vast majority of 
low-dimensional materials realize instances of the so-called 
Heisenberg model (see e.g. [26]). In fact, the demanding re- 
quirements to explore the Ising magnetism in real materials 
(e.g. precise one-dimensionality, strong Ising anisotropy, and 
weak exchange couplings matching the available magnetic 
fields) have postponed its observation until the recent exper- 
iments with CoNb206 [27-30]. We stress that prior to this 
experiment, trapped-ion QSs had already targeted the onset 
of Ising criticality [8, 10]. More recently, Ising interactions 
in a large collection of ions have been observed in Penning 
traps [12], which in combination with the recent results for 
neutral atoms [31], yield the unique possibility of tailoring 
the microscopic properties of the quantum Ising magnet (e.g. 
couplings, dimensionality, geometry). 

In this manuscript, we investigate the capabilities of trapped 
ions as QSs of frustrated quantum Ising models (FQIMs). 
Here, frustration arises from the impossibility of minimiz- 
ing simultaneously a set of competing commuting interactions 
(i.e. classical frustration) [32], whereas the quantum fluctua- 
tions are introduced by a transverse magnetic field. We remark 
that, with the exceptions of three-dimensional spin ice [33] 
and some disordered spin-glass compounds [34], most of the 
frustrated materials correspond once more to Heisenberg mag- 
nets [35]. This would turn the extensive research on low- 
dimensional frustrated Ising models [36, 37] into a pure theo- 
retical enterprise. Fortunately, the seminal experiment [9] has 
proved the contrary by demonstrating that the physics of small 
frustrated networks [17] can be explored in trapped-ion labo- 



ratories. In a recent work [18], we have studied a different 
approach to control the degree of frustration in ion crystals, 
which works independently of the number of ions and is thus 
amenable of being scaled to large systems. In this manuscript, 
we elaborate on that proposal to develop a versatile quantum 
simulator for a variety of spin ladder geometries. 

A quantum spin ladder is an array of coupled quantum 
spin chains. Each of these spin chains is usually known as 
a leg of the ladder, whereas the couplings between them de- 
fine the ladder rungs. Let us note that the study of antifer- 
romagnetic Heisenberg ladders has a long tradition inspired 
by the experiments with insulating cuprate compounds [38], 
which become high-temperature superconductors after dop- 
ing [39]. These strongly-correlated systems lead to fasci- 
nating and thoroughly- studied phenomena [40]. In contrast, 
the quantum Ising ladders have remained largely unexplored, 
probably due to the absence of materials that realize them. 
Our work discusses a proposal to fill in such gap, which is 
based on techniques of current trapped-ion experiments. 

This article is organized as follows. In Sec. II, we describe 
the properties of the ladder compounds formed by a collec- 
tion of ions confined in a radio-frequency trap. We discuss 
the experimental conditions leading to a particular vibrational 
band structure, which shall be exploited to explore the physics 
of magnetic frustration. In section III, we support the above 
conclusions by the numerical simulation of a particular case: 
the trapped-ion zigzag ladder. A discussion of the experimen- 
tal feasibility of our scheme with state-of-the-art setups is also 
presented in Sec. IV. Here, we also present initial experiments 
regarding the appropriate trap design for the QS, and discuss 
possible sources of error. In Sec. V, we list the many-body 
phenomena that can be addressed with the proposed QS, rang- 
ing from ordered phases of quantum dimer models, to disor- 
dered quantum spin liquids. In Sec. VI, we address in detail 
the phase diagram of the dipolar J\-J2 quantum Ising model. 
Finally, we present our conclusions in Sec. VII. In the Ap- 
pendixes A, B and C, we discuss technical details regarding 
both the spin-dependent dipole forces and spontaneous decay, 
the micromotion, and the thermal secular motion. We describe 
conditions under which they do not affect the proposed QS. 



II. TRAPPED-ION QUANTUM ISING LADDERS 
A. Geometry of the trapped-ion ladders 



An ensemble of N atomic ions of mass m, and charge e, 
can be trapped in a microscopic region of space by means of 
radio-frequency fields [6]. This system is described by 




where {(Q a } a=XJiZ represent the effective trapping frequencies 
(see also Sec. IV and Appendix B), and {r i -,p i }^ 1 are the po- 
sitions and momenta of the ions. Note that we use gaussian 
units and fi = 1 throughout this text. 
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With the advent of laser cooling, it has become possible to 
reduce the temperature of the ions to such an extent that they 
crystallize. This gives rise to one of the forms of condensed 
matter with the lowest attained density, ranging from clus- 
ters composed by two ions [41], to crystals of several thou- 
sands [42]. Interestingly, these crystals can be assembled se- 
quentially by increasing the number of trapped ions one by 
one, so that the aforementioned transition to the many-body 
regime [1] can be studied in the laboratory. Besides, the geom- 
etry of the crystals can be controlled experimentally by tuning 
the anisotropy of the trapping frequencies [43], which opens 
a vast amount of possibilities for the QS of magnetism. 

The ion equilibrium positions {r®}f =l are determined by 
the balance of the trapping forces and the Coulomb repulsion. 
Formally, they are given by V r .V = 0, where V contains the 
trapping and Coulomb potentials. By introducing the unit of 
length l z = (e 2 /mco^y/ 3 , the dimensionless equilibrium posi- 
r? / l z follow from the solution of the system 
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where we have introduced the anisotropy parameters K a = 
(cQ z /cQ a ) 2 - As shown in the experiments [44], by increasing 
the parameters k x = K y towards a maximum value K* max = 1, 
the geometry of the ion crystal undergoes a series of phase 
transitions starting from a linear chain, via a two-dimensional 
zigzag ladder, to a three-dimensional helix. We remark that 
these structural phase transitions have been the subject of con- 
siderable interest on their own [45, 46]. 

We have recently investigated the possibility of setting a 
large anisotropy between the radial trap frequencies 1 > K x ^> 
K y , such that the ion crystals get pinned to the xz-plane [18]. 
This requires a modification of the symmetric electrode con- 
figuration of the usual linear traps (see the details in Sec. IV). 
Besides, the parameter k x can be increased by a bias voltage, 
yielding a variety of geometries. For instance, above a certain 
value K* C5 2, the linear chain transforms onto a 2-leg ladder (see 
Figs. l(a)-(b)). By solving numerically the system of equa- 
tions (2), we observe that the 2-leg zigzag ladder first evolves 
into a 3 -leg ladder as k x is increased, then yields a 4-leg lad- 
der [see Figs. 1(c) and 1(f)], and so on. Hence, there should 
be a sequence of structural transitions at the critical values 
Kc,2 < *Q,3 < Kc,4 < • • • < K"c,rcp where n\ determines the num- 
ber of legs of the trapped-ion ladders. The rungs of these lad- 
ders correspond to the diagonal links, yielding a collection of 
bond-sharing triangular plaquettes. In Sec. IV, we describe 
a method to build ladders from corner- sharing triangles (see 
Figs. l(d)-(e)), which widens the applicability of our QS. 

Eventually, when K x —> 1, the crystal must correspond to 
a two-dimensional bond-sharing triangular lattice, or corner- 
sharing Kagome lattice, with ellipsoidal boundaries. Accord- 
ingly, not only can we control the number of legs of the 
trapped-ion ladder, but also explore the crossover from quasi- 
one-dimensional physics to the two-dimensional realm. With 
respect to the FQIMs, this two-dimensional limit is very ap- 
pealing due to its connection to quantum spin liquids [24], and 
quantum dimer models [37]. However, as happens for Heisen- 
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Figure 1. Trapped-ion ladders: Self- assembled geometries for 
N trapped ions as a function of the anisotropic trapping frequencies 
(Oy > co x , co z : (a) Linear chain for N = 50 and K x = 4 • 10 -4 , (b) Two- 
leg zigzag ladder for N = 51 and K x = 1 • 10 -2 . (c) Three-leg ladder 
of bond- sharing triangles for N = 52 and K x = 4 • 10 -2 , (d) Two-leg 
ladder of corner- sharing triangles for N = 52. From this geometry, it 
is possible to obtain the two following ladders by selectively hiding 
some of the ions (grey crosses) following the methods outlined in 
Sec. IV. (e) Three-leg ladder corresponding to a Kagome stripe of 
corner- sharing triangles and hexagons for N = 52, (f) Four-leg ladder 
of corner- sharing triangles for N = 53 and K x = 8 • 10 -2 . 



berg magnets [38], the ladder compounds may already contain 
fascinating phenomenology. Let us remark that the precise 
knowledge of the crystal plane is a fundamental advantage for 
the quantum simulation of FQIM [18]. Prior to the discussion 
of the effective FQIM, we will describe the collective vibra- 
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tional modes in these ladder compounds, since they shall act 
as mediators of the magnetic Ising-type interaction. 

At this point, it is worth commenting on the possibility of 
engineering the lattice positions of the ion crystal by means of 
micro-fabricated electrode arrangements, (e.g. surface traps 
as proposed in [47]). If the technical problems related to the 
anomalous heating observed close to the electrodes are over- 
come, these new generation of traps could be combined with 
the proposed QS to study frustrated spin models in arbitrary 
lattices without the complication of micromotion. However, 
to keep within reach of the current technology, we focus on 
the more standard Paul traps where experiments on quantum 
magnetism have already been performed [8-10]. 



B. Collective vibrational modes of the ion ladder 



Due to the Coulomb interaction, the vibrations of the ions 
around the equilibrium positions, r; = r? + Ar,, become cou- 
pled. By expanding the Hamiltonian (1) to second order in the 
displacements Ar,-, one obtains a model of coupled harmonic 
oscillators corresponding to the harmonic approximation 
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In this expression, the coupling matrix between the harmonic 
oscillators can be expressed in terms of the mutual ion dis- 



tance r?- = r? - 
of ions Q?P = 



r°j, and the quadrupole moment for each pair 
-e[3(r°.)«(if> " M r ° ) 2 ] as follows 
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For the particular ladder geometries in Fig. 1, we identify 
two types of vibrational excitations, the so-called transverse 
modes that correspond to the vibrations of the ions perpen- 
dicular to the ladder (i.e. v-axis), and the planar modes that 
account for the coupled vibrations within the ladder (i.e. jc, z 
axes). The transverse modes are decoupled from the planar 
vibrations, and are described by a set of coupled oscillators 
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where COy ->> &f = CO y {\ - KyL/^- |rf - rf | 3 ). This term 
amounts to the Einstein model of individual lattice vibrations, 
where the renormalization of the frequencies is caused by the 
mean-field-type interaction of one ion with the rest of the ion 
ensemble. In our case, we must also consider the coupling be- 
tween distant oscillators, whose magnitude relative to the trap- 
ping frequencies scales as (e 2 / 7|) / (mco 2 ) = K y <C 1. Hence, 
the transverse vibrations are described by a set of harmonic 
oscillators with weak couplings that decay with a dipolar law. 

The planar vibrations are more complex since the mo- 
tion along both axes, x and z, becomes coupled through the 



Coulomb interaction. In this case, the Hamiltonian is 
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where 6bf x = 00% (1 + KxL/^£^|if - f°| 5 ) (equivalently for 
z), and Qff is expressed in terms of the dimensionless equilib- 
rium positions. Since we have assumed that K x ^> K y , the pla- 
nar vibrations are coupled more strongly than the transverse 
ones. In the last term of this Eq. (6), we observe how the non- 
diagonal terms of the quadrupole, ^ 0, are responsible for 
the coupling of the motion along the x and z axes of the ladder. 

So far, we have reduced the original Hamiltonian (1) to a 
pair of quadratic boson models (5)-(6). These can be exactly 
solved by introducing the canonical transformations 

N 1 



2N | 

n =i y2mkl n 



(J) 



2N | 

n =\ \/2mil n 

where a\,a n stand for the N creation- annihilation operators 
for the quantized excitations (i.e. phonons) of the transverse 
modes, and b\^b n are the the corresponding 2N operators for 
the planar modes. In these expressions, deter- 
mines the amplitude of the transverse (planar) oscillations of 
the ion at site i due to the collective vibrational mode labeled 
by n, such that Q^-(Q^) are the corresponding normal-mode 
frequencies. Formally [48], j^^ n are given by the or- 

thogonal matrices that diagonalize the second order Coulomb 
couplings in Eqs. (5)-(6), or equivalently 
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where f ; f P 



;/ , l} , } are dimensionless couplings. 

These equations (8) must be solved numerically with the pre- 
vious knowledge of the equilibrium positions (2), and yield 
the following quadratic phonon Hamiltonian 
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where the index n labels the normal modes with an increasing 
vibrational frequency > > 
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Figure 2. Vibrational modes of the ladder compounds: (a) 

Schematic representation of the phonon branches, whereby the trans- 
verse vibrational frequencies £l„ span around the trap frequency COy, 
and are situated far away from the planar phonon branch Q!^, which 
spans around co x , co z . (b) Laser beam arrangement (inset) for a spin- 
dependent dipole force with an effective wavevector kL in the xy- 
plane slightly tilted from the v-axis. The effective frequency &>l is 
tuned above the resonance of the transverse modes (see (a)), such 
that the detuning s fulfill 8 U ^> 8±. 



We now discuss a qualitative picture for the phonon 
branches valid for all the different ladders. As argued above, 
the transverse phonon modes are weakly coupled due to the 
small parameter Ky<L By inspecting Eq. (8), one realizes 
that this condition will lead to a transverse phonon branch 
Q^- G [£2^,^] with a small width around the trap frequency 
COy (see Fig. 2(a)). This property is not fulfilled by the pla- 
nar modes Q}„ G , Q^] , which are more strongly coupled 
K x ^> K y , and will thus present a wider branch around (Q X ,CQ Z . 
However, due to the frequency anisotropy COy ^> CQ X ,CQ Z , the 
planar phonon branch will always be situated far away from 
the transverse-mode frequencies, even if its width is consider- 
ably larger. This property will turn out to be essential for the 
QS of frustrated spin models [18], whereby the role of the spin 
is played by two electronic levels of the ion, and the interac- 
tions are mediated by the collective vibrational excitations. 

The main idea is that such a clustered phonon spectrum will 
allow us to use a pair of laser beams that only couple the spins 
to the transverse phonon modes, even when the radiation re- 
sulting from the interference also propagates along the plane 
of the ladder (Fig. 2(b)). The interest of this idea is two-fold. 
On the one hand, the transverse phonon modes are the ideal 
mediators of the spin- spin interactions due to their higher in- 
sensitivity to ion-heating and their lower contributions to the 
thermal noise in the QS, as compared to the planar modes. 
On the other hand, by using a laser configuration with a com- 
ponent along the rungs of the ladder, we can exploit the ra- 
tio of its effective wavelength with the ion mutual distances 
in order to tailor the sign (ferromagnetic/antiferromagnetic) 
and the magnitude of the spin- spin interactions anisotropically 
(i.e. depending on the direction joining the pair of ions). As 
shown below, this opens the possibility of realizing a versatile 
QS of frustrated quantum spin ladders, which is amenable of 
being scaled to larger ion numbers and a variety of geometries. 



C. State-dependent dipole forces 

The use of the collective vibrational modes as a com- 
mon data bus to perform two-qubit gates can also be under- 
stood in terms of phonon-mediated spin-spin interactions [49], 
which opens the route towards the QS of quantum mag- 
netism [15, 16]. We now introduce the key ingredient for such 
QSs, namely, a spin-phonon coupling that originates from a 
laser-induced state-dependent dipole force [50, 51]. 

So far, our discussion applies to all ion species and radio- 
frequency traps. However, in order to exploit the phonons as 
mediators of a spin-spin interaction that is sufficiently strong, 
the trap frequencies for each particular ion species must be 
tuned such that l z = (e 2 /mco 2 ) 1 / 3 ^ 1-10 /im, and COy > 0) x > 
co z . This typically constraints the order of magnitude of the 
trap frequencies to the range co z /2n, co x /2n ^ 0.1-1 MHz, and 
£0y/27r«l-lO MHz. Moreover, since we aim at a flexible con- 
trol of the anisotropy of the spin couplings, we shall focus on 
the ion species that allow for a two-photon lambda scheme to 
implement the spin-phonon coupling. Hence, our discussion 
will be specific to singly ionized alkaline-earth ions, either 
with a hyperfine structure 9 Be + , 25 Mg + , and 43 Ca + , or with a 
pair of Zeeman-split levels such as 40 Ca + and 24 Mg + [52]. 

In Fig. 3, we represent schematically the atomic energy lev- 
els of such ions, which have a single valence electron in the 
orbital n 2 5 , 1 / 2 that can be optically excited to n 2 P\/2i n2 P?>/2 
via a dipole-allowed transition. Here, we use the standard no- 
tation n 2iS+1 L/, with n as the principal quantum number, and 

L, / as the spin, orbital, and total electronic angular momen- 
tum. Depending on the nuclear spin / of the particular ion, 
the ground- state manifold will be split into a pair of Zeeman 
states (1 = 0), or a set of hyperfine levels (7^0). We select 
two of such states \U), which are separated by an energy 
gap C0o, to form the effective spins of our QS (see the inset of 
Fig. 3). For hyperfine spins, the energy splitting is C0q/2tc w 1- 
10 GHz, whereas for Zeeman spins its order of magnitude de- 
pends on the external magnetic field cOq/Bq « 271 x 10 GHz/T. 

Accordingly, the resonance frequencies usually lie in the 
radio-frequency/microwave regime. However, it is not cus- 
tomary to use radio-frequency/microwave radiation to couple 
the spins to the collective vibrational modes, since its wave- 
length is much larger than the typical ion oscillations (but see 
the recent experimental progress [53]). A possible alterna- 
tive is to use a pair of laser beams with optical frequencies 
ft>i, 0)2, a "moving standing wave", to induce a two-photon 
stimulated Raman transition through an excited state |r t -) (see 
Fig. 3). When the laser beams are far off-resonant with respect 
to this dipole-allowed transition, such that the detuning A is 
much larger than the decay rate of the excited state T and the 
laser Rabi frequencies s,^2,s, where s =t,>U it is possible 
to eliminate the excited state from the dynamics and obtain 
a Hamiltonian that only involves the spins and the phonons 
(see Appendix A for the effective master equation). In partic- 
ular, when the laser beatnote is C0\ = C0\ Oh ^ co y <C C0q <C A 
(Fig. 2(a)), one obtains the following ion-laser Hamiltonian 
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Figure 3. Lambda scheme for the dipole force: Two electronic 
states { 1 1/) , | ij) } with an energy difference (Oq are selected from the 
ground-state manifold n 2 Si/ 2 to form an effective spin- 1/2. These 
states can be manipulated by a pair of laser beams with Rabi fre- 
quencies £2i ?C 7, ^2,a mat induce a transition to an excited state \rj) in 
the n 2 P^j2 manifold from the spin state s =t,i m n2 S\/i- Note that 
the transition wavelength A sp lies in the optical regime. When the 
beatnote co^ = C0\ — (O2 is tuned close to the trap frequency C0y, such 
that the detuning A is much larger than the spontaneous decay rate 
r of the transition and the Rabi frequencies, one obtains the desired 
state-dependent dipole force. 



which can be interpreted as a time-dependent differential ac- 
Stark shift due to the interference of the two laser beams. 
Here, we have introduced of = |ti)(t/| _ ll/)U/|> an d the 
two-photon differential Rabi frequency £2l = ^2 1 ~ 
£^ ^£2|^)/2A. The corrections due to the spontaneous decay 
from the excited level contribute to the Rabi frequency with 
^ L + (r/A) 2 ), and can be thus neglected if T <C A. 

The idea is that by controlling experimentally the polariza- 
tions, intensities and detunings of the laser beams, one can 
find regimes where the effective Rabi frequency £2l does not 
vanish and leads to a spin-dependent dipole force. Let us note 
that the adiabatic elimination also leads to standard ac- Stark 
shifts and to a running wave that only couples to the vibra- 
tional excitations. These terms must be compensated by care- 
fully selecting the laser-beam parameters. Let us finally em- 
phasize that the effective Raman wavevector kL = ki — k2 = 
&l cos 6e x + &l sin 6e y has a component along the ladder plane 
(see Fig. 2(b)), and in principle couples to both the planar 
and transverse phonons. The reason for the announced decou- 
pling from the planar vibrational modes is that the beatnote of 
the laser beams, which is near the resonance of the transverse 
phonons, will lie far off-resonance with respect to the planar 
vibrational modes (Fig. 2(a)). This qualitative argument will 
be quantified below, and supported numerically in Sec III. 

After introducing the phonon operators in Eq. (7), we per- 
form a Taylor expansion of Eq. (10) for the small transverse 
and planar Lamb-Dicke parameters 
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By setting COl > CQy, such that the bare detuning fulfills \S y \ = 
\cOy — COl I <^ \8 X \ = \co x — col |, ®y, one can neglect all non- 
resonant terms apart from a state-dependent dipole force that 
couples the spins to the transverse phonons. In the interaction 
picture with respect to the phonon Hamiltonian (9), we get 
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where S n ± — £l„ — COl- We note that, in order to neglect all 
the remaining terms of the Taylor expansion, a rotating wave 
approximation (RWA) must be performed, provided that 



col 



«1, 



(13) 



The first condition is required to neglect the off-resonant con- 
tributions to the ac- Stark shift of the energy levels. Addi- 
tionally, by virtue of Eq. (11), this condition ensures that the 
spin-phonon couplings involving higher-order powers of the 
transverse phonon operators can also be neglected. For the 
regime considered in this work, namely COl ~ C0y ^> co x > CO z , 
it will suffice to consider Q^jln < 0.1-1 MHz to accomplish 
this constraint. The second condition in (13) is necessary to 
avoid that the dipole force also couples the spins to the planar 
phonons. Hence, it characterizes the parameter regime where 
our previous qualitative discussion about the decoupling of 
the planar vibrational modes holds. The fulfillment of this 
condition relies on the large gap between the frequencies of 
the transverse and planar normal modes. Besides, by setting 
the laser-beam arrangement such that 6 w n/2, one obtains 
r\ m <C 1, which warrants the fulfillment of the last condition. 
In Sec. Ill, we will confirm the validity of these constraints 
numerically for the particular case of a zigzag ion ladder. 

Depending on the spin state, the dipole force in Eq. (12) 
pushes the ions in opposite directions transversally to the lad- 
der plane. Besides, the phase of this pushing force depends 
on the ratio of the ion equilibrium positions and the effective 
wavelength of the interfering laser beams. As shown below, 
this is precisely the parameter that shall allow us to control the 
anisotropy of the effective spin- spin couplings. 



D. Spin models with tunable anisotropy and frustration 

There are numerous situations in nature where interactions 
are mediated by the exchange of particles. Of particular rele- 
vance to the field of magnetism is the so-called Ruderman- 
Kittel-Kasuya-Yosida (RRKY) mechanism [54], whereby a 
Heisenberg coupling between distant nuclear spins is medi- 
ated by electrons from the conduction band of a metal. The 
sign of the Heisenberg couplings alternates between ferro- 
magnetic/antiferromagnetic as a function of the ratio between 
the Fermi wavelength and the mutual spin distance. We have 
recently shown [18] that a similar phenomenon occurs for 
trapped-ion ladders subjected to the dipole force in Eq. (12). 
Instead of a Heisenberg-type interaction between the nuclear 
spins, one obtains a periodically-modulated Ising-type cou- 
pling between the spins formed by two electronic states of the 
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ion. Interestingly enough, the modulation can be experimen- 
tally tailored by controlling the direction of propagation of the 
interfering laser beams providing the dipole force. 

The Ising interaction between the effective spins of two dis- 
tant ions can be understood as a consequence of the virtual 
phonon exchange between these ions. The dipole force (12) 
pushes the ions transversally, exciting thus the transverse vi- 
brational modes. These phonon excitations, being collec- 
tive, can be reabsorbed elsewhere in the ion crystal, provid- 
ing a mechanism to couple the distant spins. The exact ex- 
pression can be obtained by a Lang-Firsov-type transforma- 
tion [15, 55-59] that decouples the spins from the phonons 



Us- 



^Vie ikLr ?-^ 



H.c. 



(14) 



This canonical transformation leads to an effective Hamilto- 
nian Ul (B p + H d )U s « tfeff + Bp, where we consider a pic- 
ture such that the vibrational Hamiltonian H v (9) absorbs the 
time-dependence of the dipole force (12), namely | \jf(t)) = 
U(t)\ y(*))> wnere U(t)= exp(k£ n co^ala n ). This leads to the 
following representation of the total Hamiltonian 

B p +B d = £5„ x ak + ^Z^'Vn^itofat + H.c). 

(15) 

After applying the above Lang-Firsov-type transformation, 
and moving back to the original Schrodinger picture, we ob- 
tain the following effective Ising model 



(16) 



The effective spin couplings have the following expression 



zeff 



y ft^sin 2 



cos(k L -r°,) (17) 



As announced previously, the phase dependence of the dipole 
force (12) on the ratio of the ion equilibrium positions and 
the effective wavelength of the light has been translated in 
the particular periodic modulation of the interaction strengths 
JfJ* oc cos(kL • r? ). This sign alternation is similar to that 
found in RKKY metals, such that the transverse phonons play 
the role of the conduction electrons, and the wavelength of the 
interfering laser beams acts as the Fermi wavelength. In the 
regime of interest for the ladders, jc y <l, the couplings 



ff J ef f COS <j>ij 
J ij | ? 0_ ? 0|3^ 



(18) 



display a dipolar decay law, where we have introduced 

^eff=^#W fo = 2ft / , (19) 



852 



such that T] y = &l sin 6 / ^Jlmddy is the bare Lamb-Dicke pa- 
rameter. From these expressions, it becomes apparent that the 
interactions between spins belonging to the same leg of the 



ion ladder (fyj = 0) correspond to antiferromagnetic JfJ f > 
Ising couplings. Conversely, the interactions between the 
spins from different legs of the ladder (fyj ^ 0) can be fer- 



romagnetic Jfj < or antiferromagnetic Jfj* > depending 
on the laser parameters. We can thus tune the sign and mag- 
nitude of the spin-spin couplings anisotropically. Note that, 
even if the typical optical wavelengths are much smaller than 
the mutual ion distances Al <C l z « 1-10 fim, the angle 6 can 
be tuned around 6 « | so that fyj attains any desired value 
<j>ij G [0, 2n] . Alternatively, it is also possible to maintain the 
laser-beam arrangement fixed, and modify the anisotropy ra- 
tio K x = (co z /co x ) 2 in order to control the inter-ion distances 

— x°j, attaining thus the desired value fyj G [0,2;r]. Either of 
these two methods will turn out to be essential to explore the 
full phase diagram of the frustrated quantum spin ladders. 

This anisotropy of the model leads to frustration when 



-J) 



eff I 



(20) 



where P stands for the elementary plaquette of the lattice, 
which is a triangle in our case = — 1 (see Figs. 1(b), (c) 
and (e)). Therefore, the Ising model is frustrated if there is an 
odd number of antiferromagnetic couplings per unit cell. We 
note that this standard criterion of geometric frustration [60] 
can be extended to situations where quantum fluctuations and 
frustration have the same source [61]. In our case, however, 
the frustration is purely classical, interpolating between the 
antiferromagnetic frustration, and the frustration due to com- 
peting ferromagnetic and antiferromagnetic interactions. In 
order to introduce quantum fluctuations, a microwave directly 
coupled to the atomic transition yields 

H^^^Jffafa^-h^af, (21) 

i j^i i 

where we have introduced of = (|;| + \U) (t/|» an d h plays 
the role of an effective transverse field of strength h due to the 
microwave, which is responsible for the quantum fluctuations. 

Let us close this section by rewriting the effective 
anisotropic quantum Ising model (21) in a notation that is 
more appropriate for quantum spin ladders (see Fig. 4). We 
substitute the label of the spins i = 1, • • • ,N for two new in- 
dices that account for the number of legs y = 1, • • • ,n\, and 
the number of spins in each of these legs i s = 1 , • • • , Ly, such 
that £ r L r = N. Then, the Hamiltonian can be rewritten as 
// e ff = //leg + H mng , where the leg and rung Hamiltonians are 



(22) 



Here, we have introduced the intra- and inter-leg Ising cou- 
pling strengths, which are respectively 

•Jeff 



Jett 



(23) 
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Figure 4. Anisotropic spin-spin interactions: The legs of the 
trapped-ion ladder are labeled by the index y = 1 , . . . , n\, whereas the 
spins in each leg correspond to i s = 1 , . . . , Ly. The effective spin lad- 
der hamiltonian is composed of two terms: H\ eg contains single- spin 
terms corresponding to the transverse field h ^ (yellow-red laces) 
and the antiferromagnetic dipolar couplings J > (yellow bonds). 
// m ng contains the couplings between spins of different rungs, such 
that / ^ (yellow-red bonds) depends on the mutual rung distance 
dyfx. We also show the unit vectors ai , a2 of the triangular lattice. 



III. A DETAILED CASE: THE ZIGZAG LADDER 

A. Numerical support for the anisotropic Ising model 

We consider the simplest scenario where the anisotropy of 
the Ising model can be tested, namely, a three-ion chain in 
a zigzag configuration. We consider the following guiding 
numbers for the trap frequencies (Qy/(Q Z = 20, (Q x /co z = 1.43, 
and (O z /2n = 1 MHz, although we emphasize that the scheme 
will equally work for different values as far as the above con- 
straints are met. These values lead to the equilibrium positions 

rf e {(-0.22,0, -0.92), (0.44,0,0), (-0.22,0,0.92)}. (24) 

We focus on the crucial assumption that allows us to derive 
the effective frustrated Hamiltonian (21), namely the possi- 
bility to neglect the pushing force on the planar modes. To 
address the validity of this approximation, we start from the 
vibrational Hamiltonian in Eq. (3), and introduce the creation- 
annihilation operators for the local ion vibrations 



where 7 e ff = / e ff cos (j) r{1 , such that (j) r!1 = &l cos 0d r{1 , and d r{1 
is the inter-leg distance (see Fig. 4). These are the central 
equations of this manuscript, describing a general n\-leg quan- 
tum Ising ladder, whereby the each of the legs corresponds 
to an antiferromagnetic quantum Ising chain with long-range 
dipolar interactions. The one-dimensional chains are coupled 
by means of dipolar Ising pairwise interactions with a strength 
and sign that can be experimentally controlled. In particular, 
we shall be interested in a ferromagnetic coupling that com- 
petes with the intra-leg antiferromagnetic interactions, leading 
thus to the phenomenon of magnetic frustration. 

At this point, it is worth commenting on the interesting re- 
cent proposal [62] for the simulation of any network of N in- 
teracting spins by using a linear crystal of ions. By exploiting 
N different laser beams individually addressed to each ion, 
such that the detunings and Rabi frequencies are fine tuned, it 
is possible to synthesize the connectivity of any desired net- 
work. Our approach is different since it exploits the specific 
geometry of self-organized planar ion crystals, and can be 
scaled to large ion ensembles straightforwardly. Additionally, 
it benefits from the simplicity of using a single laser-induced 
dipole force that lies far off-resonance from the whole vi- 
brational branch. This contrasts the proposal in [62], where 
the forces lie within the vibrational branch, and resonance ef- 
fects must be carefully avoided for larger ion chains where 
the phonon branches become denser. We note that the scheme 
in [62] has a higher flexibility in the simulated lattices, al- 
though the geometry of the ladders in our approach can be 
partially modified following the prescriptions of Sec. IV. 

In the following section, we support the validity of this an- 
alytical treatment with a numerical study of the zigzag ladder. 



(25) 

Even if the state-dependent dipole coupling in Eq. (10) only 
acts along the xy plane, the motion along the z-axis gets cou- 
pled through the Coulomb interaction (6). Therefore, we must 
treat the complete vibrational Hamiltonian 



i a=x,y,z 
CO 



(26) 



ij 



The state-dependent dipole force can also be written in this 
local basis, yielding the following Hamiltonian 

#2=i^E E ^haOfale-^'+H.c, (27) 

^ i a=x,y 



where r\ a = k\^ a / ^2m(O a are the bare Lamb-Dicke factors. 

In order to integrate numerically the Schrodinger equation 
for the timescales of interest tf w l/|/ e ff|> which are three or- 
ders of magnitude larger than the timescale set by the dynam- 
ics of the dipole force 1/cOl, it would be desirable to work in 
a picture that absorbs the fast time-dependence of (27). This 
is possible if we neglect the counter-rotating terms of Eq. (26) 
that correspond to phonon non-conserving processes, which 
is justified by a rotating- wave approximation when 



co- 



(28) 



Then, it is possible to move to a picture where the 
creation-annihilation operators rotate with the laser frequency, 
\\j/(t)} = exp(k^ a; 0)L^J a (2- a )|i/A(^)), and the Hamiltonian 
that will be numerically explored H = H\ +H2 becomes time- 



independent, namely 

fix =E 5 ««U« +?E E(^) 1/4 ^f 4« y/3 , 



2 



a,j3 i'j 



(29) 



Note that this unitary transformation does not affect the spin 
dynamics, and is thus well- suited to study the validity of our 
previous derivation of the effective Ising model. 

In order to account for two sources of noise that are usually 
the experimental limiting factors for spin-oriented QSs, we 
include a fluctuation of the atomic resonance frequency 



5 3 =£§A£(f)of, 



(30) 



where Ae(t) is a stochastic process. This term may corre- 
spond to the Zeeman shift of non- shielded fluctuating mag- 
netic fields, or to a non-compensated ac-Stark shift caused by 
fluctuating laser intensities. Its dynamics can be modeled as a 
stationary, Markovian, and Gaussian process [63], as follows 



Ae(t + St) = Ae(f)e _ * + [f (1-e" 



(31) 



where n g is a unit Gaussian random variable, and c, T char- 
acterize the diffusion constant and the correlation time of 
the noise. For short correlation times T <C t, one obtains an 
exponential damping of the coherences with a typical time 
72 = 2/ct 2 . We set r = 0.172, and 72 « 10 ms, which is a 
reasonable estimate for the observations in experiments. Note 
that this dephasing timescale still allows for the observation 
of the faster coherent spin dynamics l/Jff « 1 ms <C 7^. 

We study numerically the time evolution under the total 
Hamiltonian H = £ m // m in Eqs. (29)-(30), and compare it to 
the effective description 77 e ff for a vanishing transverse-field 
in Eq. (21). We note that the following numerical simulations 
focus on a ground-state-cooled crystal, where the nine vibra- 
tional modes have one excitation at most. To consider the 
dephasing noise, we integrate over N = 10 3 different histories 
of the fluctuating frequency (31), and perform the statistical 
average of the dynamics. The effects of finite temperatures, 
together with other error sources, are addressed in Sec. IV. To 
observe a neat hallmark of the anisotropy due to the modula- 
tion of the interaction strengths JfJ f °< cos(kL • r?-), we study 
the dynamics of the spin state \\jf s ) = |+)i (g> |— )2 <S> \—)3, 
where |±) = (| t) ± 1 1))/ y/2, for two sets of parameters. 

( i) Allowed spin hopping: We set the following parameters 
ri y = 0.1 = 10t] x , col = l.lctfy, and £} L = 0.15\8y\/r]y, and di- 
rect the laser beams so that 6 = n/2 (i.e. e x kL = 0). In 
this case, the interfering radiation does not propagate along 
the triangle plane, and there is no modulation of the sign of 
the couplings. We observe that the initial spin excitation |+) 
located at site 1, can tunnel to the two remaining sites of the 
triangular plaquette as a consequence of the Ising-like cou- 
pling. In Fig. 5(a), we compare the numerical results with the 
effective description. From this figure, we can conclude that 
the effective Ising Hamiltonian yields an accurate description 
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Figure 5. Spin dynamics in the triangular plaquette: (a) Dynam- 
ics of an initial spin excitation |+)i in the allowed-hopping regime. 
The dynamics of (&j(t)) obtained from the effective description H Q ff 
((of (f)) red line, (of (t)) blue line, (of (t)) yellow line) show a re- 
markable agreement with the numerical simulation of the complete 
Hamiltonian H ({of(t)) red circles, {0% (t)) blue triangles, (o$(t)) 
yellow squares), and both describe the transport of the spin excitation 
around the plaquette. (b) In the inhibited hopping regime, one also 
observes a complete agreement between both descriptions, which de- 
scribe how the spin excitation cannot occupy site 2 as a consequence 
of the vanishing Ising interaction. 



of the spin dynamics in a timescale tf ~ l// e ff~ 2// e ff ~ 1-2 
ms< 72 « 10 ms. For longer timescales, the dephasing leads 
to a larger deviation from the effective description. Note that 
this result supports the validity of the isotropic Ising interac- 
tion, but we still have to address whether our scheme to con- 
trol interaction anisotropy is also accurate. 

(ii) Inhibited spin hopping: We use the same parameters 
as before, but now consider that the interfering laser radiation 
also propagates along the plane defined by the triangle. In 
particular, we set the angle and the effective wavelength as 



follows A-l/ cos0 = 4/ z (iy —^2) = 4/ z (^3 — j^), which leads 
to the factors 0i2 = (^32 = tt/2, such that the effective Ising 
couplings between sites 1-2 and 2-3 completely vanish Jf^ = 
/|f = 0. Therefore, we reach a highly anisotropic situation 
where the spin excitation can only hop between sites 1^3. 
This is confirmed by the numerical results in Fig. 5(b), where 
a good agreement with the effective description is displayed 
once more in the timescale tf& 1-2 ms. 

Let us finally stress that an experiment with three ions 
in this triangule would be a neat proof-of-principle to show 
that anisotropic Ising models can be realized following our 
scheme. We also stress that this scheme is amenable of being 
scaled to larger systems, since the periodically modulated cou- 
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Figure 6. Spin frustration in the zig-zag ladder: (a) Scheme of 
the Ising interactions between neighboring ions in the zig-zag ladder, 
whereby the red links stand for a ferromagnetic coupling, whereas 
the yellow-red links can be antiferromagnetic/ferromagnetic inter- 
actions, (b) Frustration in a triangular plaquette where one of the 
bonds cannot be satisfied by any spin configuration (left: purely anti- 
ferromagnet, right: competing interactions), (c) Effective interaction 
strength between the central ion jo = N/2 and its neighbors, as ob- 
tained from the numerical solution of Eq. (17) for N = 30 ions in a 
trap with C0y ^> co x = 6.1co z . For 0/ o j o +i = 0, one observes the dipo- 
lar decay of the frustrated antiferromagnetic interactions, (d) For 
0/o ,70+ 1 =7i/2, the sign alternation leads to competing interactions. 



plings are controlled globally by the ratio of the laser wave- 
length to the inter-leg equilibrium positions, and thus does not 
depend critically on the size of the Coulomb crystal. 



B. Frustration by competing interactions 



Once the validity of the effective Ising model (21) has been 
numerically supported, we can exploit the anisotropic inter- 
actions (17) to interpolate between a frustrated Ising ladder 
due to antiferromagnetic couplings, or due to the competi- 
tion of ferromagnetic and antiferromagnetic interactions. In 
Fig. 6(a), we present a scheme of the spin frustration. The 
intra-leg spin couplings in the ladder (yellow) correspond to 
antiferromagnetic interactions which, according to Eq. (17), 
cannot be modulated. On the other hand, the inter-leg cou- 
plings along the diagonal rungs of the ladder may correspond 
to antiferromagnetic (yellow) or ferromagnetic (red) interac- 
tions depending on the value of 0; 7 - = &lcos 0(*9 — Xj). Both 
situations lead to frustration (see Fig. 6(b)), since only two of 
the bonds can be satisfied simultaneously. By using the nor- 
mal modes of the inhomogeneous zigzag chain, we compute 
numerically the spin couplings, and show that for <t>j j +i = 0> 
we obtain an antiferromagnetic coupling (Fig. 6(c)), whereas 
an alternating sign arises for 0/ o j o +i = 71 / 2 (Fig. 6(d)). 



IV. EXPERIMENTAL CONSIDERATIONS 

In previous sections, we discussed the regimes of validity 
of the spin models by both analytic and numerical methods. 
Here, we analyze the capability of ion-trap experiments to 
meet the required conditions, considering current technology 
limitations and possible sources of error in ion-trap experi- 
ments. This study is supported by initial experiments. 



A. Trap design study 

A crucial condition for the frustrated quantum spin mo- 
dels is COy ^> co x > co z , which leads to the clustering of vi- 
brational branches exploited in Sec. II. This condition implies 
that the rotational symmetry of the trap potential in the xy- 
plane, which is a common property of linear Paul traps, must 
be explicitly broken. In this section, we discuss two possi- 
ble strategies to achieve this goal, and present supporting evi- 
dence based on numerical and experimental results. 

The first approach can be implemented in any linear Paul 
trap, such as the symmetric electrode arrangement shown in 
the inset of Fig. 7(a). By applying a positive offset volt- 
age £/ ffset on the DC electrodes, the trapping potential be- 
comes stronger along the diagonal direction joining the DC 
electrodes, and the trapping frequencies fulfill COy > co x (and 
vice versa for a negative offset voltage). For the experimen- 
tal results shown in Fig. 7(a), resonant radio-frequency radi- 
ation was used to excite the different modes, such that the 
induced ion motion was observed on a CCD camera. This 
method allows for the estimation of the trap frequencies, 
which show a clear agreement with numerical ion-trajectory 
simulations [64] (main panel of Fig. 7(a)). Note that the 
anisotropy of the radial frequencies (Q x /cQy can be tailored by 
the offset voltage, which also rotates the trap axes, and effec- 
tively changes the angle 6 of the laser wavevector kL with the 
x-axis (see Fig. 7(a)). This pinpoints the possibility of shap- 
ing the spin frustration according to Eqs. (18)-(19) by modify- 
ing the electrode voltages, avoiding thus the more demanding 
modification of the laser-beam arrangement. In Fig. 8, we 
show the measured positions for a. N = 11 ion crystal. We 
show how the trap-frequency anisotropy can be exploited to 
synthesize a particular 3 -leg ladder whose equilibrium posi- 
tions match perfectly the numerical predictions. 

The second strategy is to exploit a trap design with a non- 
quadrangular arrangement of the electrodes. This breaks di- 
rectly the rotational symmetry of the trapping potential (see 
the inset of Fig. 7(b)). Both the experimental data, as mea- 
sured by laser spectroscopy, and the numerics indicate a siz- 
able splitting of the radial frequencies. Note that the depen- 
dence of the angle 6 on the offset voltage depends strongly 
on the trap geometry. Unfortunately, in the present case, the 
region with the largest tunability of still coincides with the 
lowest radial anisotropies. 

In order to optimize both effects, we have designed a new 
trap (Fig. 7(c)) that meets with the special requirements of the 
QS. While the three-ion case-study may be realized with any 
state-of-the-art experimental setup, the larger scale QS with 
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planar spin systems with about 50 ions, as shown in Fig. 1, 
will require a special trap design optimized according to the 
following conditions: (a) the trapping potential should be 
highly anisotropic C0y ^> co x , (b) the trap should allow for the 
confinement of a large number of ions N < 100, (c) trap fre- 
quencies should be high such that an initialization of the ion 
crystal in a low thermal motional state is possible, (d) the trap 
geometry should reduce excess micromotion, and (e) optical 
access should allow for readout of the spin state. We now 
discuss the methods to achieve these requirements. 

The geometry for such a trap device is sketched in the inset 
of Fig. 7(c), where four round bars (r = 0.625mm) form the 
radial trapping potential. As we choose a non-quadrangular 
arrangement with distances (d Xj d y ) = (0.42 mm, 2.90 mm), 
the two simulated radial frequencies directly become non- 
degenerate. With a trapping drive frequency of £l r f/27l = 
22 MHz, and amplitude of 2.5 kV, we find a sufficiently high 
anisotropy co x /2n = 0.43 MHz< C0y/2n = 2.01 MHz along 
the requirement (a). The axial confinement is generated by 
two endcaps seperated by a distance of 25 mm, which lead 
to co z /27l = 0.23 MHz with 2 kV applied on the endcaps. In 
this potential, aiV = 4 ion crystal undergoes the first struc- 
tural transition to the zigzag configuration, and for N > 9 the 
second structural transition is obtained. We have calculated 
the equilibrium positions and vibrational modes of a three- 
legged ladder of N = 19 ions for these particular trap frequen- 
cies (Fig. 9). We note that by lowering the axial DC voltage, 
we can increase N for different crystal structures, fulfilling the 
requirement (b). Let us now address condition (c). By setting 
the trap frequency to C0y/2n =2.01 MHz, the mean phonon 
numbers after Doppler cooling lie below n y < 5. In order to 
have a thermal error below 1% for these phonon numbers (see 
Eq. (C7)), the spin-phonon coupling should be smaller than 
£2l = . 02 1 5y | / r/^ . However, this reduces the spin- spin in- 
teractions and magnetic-field noise may affect the dynamics 
at the corresponding long timescales. Therefore, multi-mode 
EIT cooling [65] to n y = 0.1 shall allow us to keep the er- 
ror rate to 1%, while maintaining spin-spin interactions in 
the kHz-regime (i.e. £2l = 0.15|5y|/rj y ). Let us stress that 
the thermal error may be minimized by considering evolution 
times that are multiples of the detuning of the closest vibra- 
tional mode. In order to estimate the effects of the excess mi- 
cromotion according to the point (d), we have estimated the 
%i values to be ^ = and \^u\ < 0.05;r with the same beam 
angles as in Eq. (B10). This should lead to relative errors on 
the order of 4-5%. Finally, the inter-ion distances are about 
10 fim for the central ions, such that we can meet with re- 
quirement (e) via high numerical aperture optics allowing for 
single-site readout. 

With the above parameters, the relative orientation of the 
ion crystal and the laser wavevector kL yields an angle of 
6 = 7r/2. As discussed above, we can vary 6 and tune the 
spin frustration by applying an offset to the DC electrodes, 
which rotates the crystal with respect to the fixed wavevector 
kL. As shown in the simulations presented in Fig. 7(c), where 
the new trap design allows for a smooth tunability of 0, while 
preserving a strong anisotropy (Q x /(dy. With £/ ffset = — 4.9 V, 
we find co x /cOy = 0.21, while the angle becomes 6 = 0.497T, 
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Figure 7. Experimental determination and numerical simulations 
of the secular frequencies in different trap geometries: Measured 
radial frequencies ca x ,C0y (green dots), calculations (green dashed 
lines), measured axial frequency co z (red dots) and calculated values 
(red dashed lines) for three different trap geometries. Furthermore 
the angle 6 (blue dots) is shown. The insets in (a)-(c) show the trap 
electrodes and axis seen from a radial cut through the trap. All exper- 
imental parameters match those of the simulation. Due to technical 
difficulties in determining the exact RF voltage we adapted its value 
by < 10%. (a) Quadrangular alignment (d x = d y = 0.82 mm) of the 
four electrodes, which leads to degenerate radial frequencies for no 
offset- voltage and only two possible values for 6. (b) Microtrap- 
design as in [66] with d x = 125jUm, and d y = 400 j um. (c) Simulated 
values for the new trap study, as described in the text. 
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Figure 8. Three-leg ion ladder: (a) Fluorescence of a 17 ion-crystal 
in the trap with (O x ^ z /{2n) = {260, 390, 111} kHz depicted in Fig. 7 
(a), inset, imaged on a CCD camera, (b) Comparison between calcu- 
lated (red cross) and measured (black circle) positions show realtive 
errors well below one percent. 



already leading to a change of sign of the spin-spin coupling 
strength Jff in Eq. (19), if we assume Al = 400 nm. 
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Figure 9. Vibrational modes and equilibrium positions of a three- 
legged ladder: (a) Equilibrium positions for a Coulomb crystal of 
N = 19 ions in a Paul trap with the axial frequency co z /2n = 229 kHz 
and the radial frequencies (O x /2n = 434 kHz,C0y/27z; = 2011 kHz. 
The frequencies are generated by a special trap (see Fig. 7(c))design, 
presented in the text, (b) Vibrational modes for the same N = 19 
ion-crystal. The transverse vibrational frequencies are clearly 
seperated from the planar vibrational frequencies Note the high 
accord with the general description in Fig. 2. 



B. Tailoring the ladder geometry 

In Figs. 1(b), (c) and (f), the ions self-organize naturally 
in a geometry of bond-sharing triangles. By controlling the 
number of legs in such triangular ladders, the QS is already 
capable of exploring a variety of interesting cooperative phe- 
nomena (see Sees. V and VI). However, it would be highly de- 
sirable to have a method to modify these geometries, widening 
thus the applicability of our quantum simulator. 

Since the geometry of the spin model (21) is determined 
by the couplings JfJ f between the spins a\ o a z - , a possibil- 
ity to tailor the geometry is to switch on/off some of these 
interactions. A possible route that allows for such a selec- 
tive coupling is a type of laser-beam hideout. Thanks to the 
relatively large distances between the ions l z ^1-10 /im, it 
is possible to address them individually with laser light [98] 
(note that the tools used so far only require global addressing). 
The idea is that, after the global optical pumping to the state 
IVs) = ®/| U)> focused laser radiation will selectively transfer 
the population to a different level \hi) that is not coupled to 
the spin-dependent force. Therefore, the ion gets effectively 
hidden. This may be achieved by polarization selection rules, 
or alternatively, by the ac- Stark shift of highly-detuned laser 
beams. Using this idea, it is possible to construct a simple 
ladder of corner-sharing triangles (see Fig. 1(c)), or a stripe 
of the two-dimensional Kagome lattice (see Fig. 1(d)). Addi- 
tionally, it opens the possibility of introducing defects in the 
lattice in order to study the role of disorder. As discussed in 
Sec. V, this tool opens many possibilities for our QS. 



C. Imperfections and noise in the quantum simulator 

Let us now address the possible sources of error in the QS. 
We place a special emphasis on the ion micromotion, which 
has not been discussed in previous QSs [8-10], since it can 
be cancelled for linear chains. We also discuss other error 
sources shared with the linear-chain QS, such as thermal fluc- 
tuations and heating of the phonons, dephasing of the spins, 
and photon scattering. We stress that, provided that the fol- 
lowing constraints to minimize the micromotion errors are ful- 
filled, the proposed ladder QS should not present more limita- 
tions than its linear-chain counterparts [8-10]. 

( i) Sources of dephasing: Usually, the most important im- 
portant source of error in the experiments is the dephasing 
of the electronic states. These might be caused by fluctuat- 
ing Zeeman shifts on the magnetic-field sensitive states re- 
quired to implement the spin-dependent dipole force (10); or 
by fluctuations in the laser-beam intensities leading to uncom- 
pensated ac-Stark shifts. A reasonable estimate of the typical 
coherence times T% w 5-10 ms [6] shows that these terms are 
important error sources in the timescales of interest tf °c 1 // eff , 
where J^/liz ~ 1 kHz. As shown in Fig. 5, the desired co- 
herent dynamics described by the effective Ising model (21) 
dominates the behavior of the system for tf « 1-2 ms. This 
sets the timescale for the QS of frustrated magnetism to be in 
the millisecond range, after which the dephasing is so strong 
that the quantum simulator is no longer faithful. 

An interesting protocol to overcome both sources of de- 
phasing simultaneously is based on the concept of continuous 
dynamical decoupling already applied to trapped ions [73]. In 
our case, it could be implemented by a strong driving of the 
carrier transition [74]. Hence, instead of relying on the spin- 
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dependent force (10), it suffices to combine a red- sideband 
term with the strong driving of the carrier. In order to obtain a 
quantum Ising model, one should introduce an additional mi- 
crowave that provides a Zeeman shift oscillating at the Rabi 
frequency of the strong carrier driving. The phase diagram can 
be explored adiabatically by an intermediate spin-echo pulse 
that refocuses the fast oscillations due to the strong driving. 

( ii) Role of the ion micromotion: The initial Hamiltonian 
in Eq. (1), upon which the derivation of the effective quan- 
tum spin ladders (22) has been built, is based on the so-called 
pseudo-potential approximation [6]. This approximation ne- 
glects the effects of the ion micromotion: a fast motion of 
frequency £1^ that is synchronous to the radio-frequency (r.f.) 
field of the Paul trap. This allowed us to focus directly on 
the slow secular motion of the trapped ions, which is de- 
scribed by an effective harmonic potential (1) with frequen- 
cies co a <C £2 r f • The validity of this assumption is well justified 
for single trapped ions, where there are experimental methods 
to minimize the effects of the micromotion [6]. For linear 
ion chains, the micromotion is also minimized when the equi- 
librium positions lie along the trap axis. Conversely, for the 
ion ladders considered in this work, the equilibrium positions 
necessarily lie off the trapping axis, leading to an additional 
micromotion that cannot be compensated. 

In Appendix B, we present a detailed discussion of the con- 
ditions under which this micromotion does not affect the QS. 
Let us comment on the possible sources of imperfections, and 
the conditions to minimize them. The r.f. heating of the trans- 
verse phonon modes responsible of the spin couplings is neg- 
ligible when 

Ky 1/2 |^fl«^, (32) 

a condition easily verified for the parameters in this work. In 
addition, during the cooling stage of the transverse vibrational 
modes, the laser frequency must be carefully tuned to avoid 
heating caused by the micromotion-induced broadening of the 
transition, or to additional micromotion sidebands [72]. 

We have also considered how the micromotion may affect 
the spin-dependent dipole force (10) and the spin-phonon cou- 
pling (12). To neglect the associated contributions, it is neces- 
sary to consider that the laser beams lie almost parallel to the 
y-axis, and are far off-resonant (see Fig. 3), namely co^/2k « 
C0y/27r «10 MHz<C a rf /27l « 0.1 GHz<C A/27T « 10 GHz. 
Besides, the Rabi frequencies of the laser beams should fulfill 

I £2 £2* I 

|ft Z|t |, « A, ht 2A U « col « £2 rf . (33) 

Finally, some care must also be placed in order not to induce 
two-photon transitions to different states \ai) of the atomic 
ground-state manifold. To avoid these processes, the relative 
Zeeman shifts of such transitions 5 fl?lS , and the associated Rabi 
frequencies must be controlled such that 

Ifti^^lS^-fflLiftrfl. (34) 

We thus conclude that, provided that the above restrictions on 
the parameters are fulfilled, the unavoidable excess micromo- 
tion of the ions in a ladder geometry should not affect our QS. 



(Hi) Thermal motion and heating of the ions: In Sec. Ill, 
we have supported the validity of the QS based on numeri- 
cal results for ground-state cooled ion crystals. However, ever 
since the early schemes for phonon-mediated gates, the ther- 
mal fluctuations of the ions have been identified as a potential 
source of errors that must be carefully considered. 

In Appendix C, we estimate the thermal error of our QS by 
both analytic and numerical methods. We solve exactly the 
Heisenberg equations of motion for the spin-dependent dipole 
force (15), and derive a scaling law for the thermal contri- 
bution to the T = spin dynamics, which is then supported 
numerically. We find that the relative thermal error, defined as 
£ t = | (of )r - (of )r=o | /| (of )r=o|, has an upper bound 



£ T < £ th = 4 — (35) 



where n y represents the mean number of transverse phonons 
in the center-of-mass mode. By considering the parameters 
used in section III, the thermal fluctuations have a small con- 
tribution (1%) for n y w 0.1 (see Fig. 19), which shows that 
perfect ground-state cooling is not required to implement the 
QS accurately. At this point, we should mention that the above 
arguments are valid for a vanishing transverse field h = 0. We 
note, however, that the contribution to the thermal error due to 
h > in the regime of interest h « / e ff has been shown [15] to 
be negligible with respect to the estimate in Eq. (35). 

In Appendix C, we also present a phenomenological master 
equation that models possible heating mechanisms in the ion 
trap. We show that the relative heating error, defined as £h = 
l(°f >ih - ( G i )r h =o|/| (of )r h =o|, where T h is the heating rate, 
only yields a small contribution (1%) to the overall error of 
the QS for heating times above 5 ms/ phonon (see Fig. 20). 

( iv) Spontaneous photon scattering: The spontaneous emis- 
sion can severely limit the advantage of quantum-information 
protocols [75-77]. In our case, the use of two electronic 
ground- states as the effective spins (see Fig. 3) makes the di- 
rect spontaneous emission | t/) ^ I U) negligible. However, 
due to the Lambda-scheme responsible for the dipole force, 
there can be photon scattering events from the intermediate 
excited state |r,-), which must be carefully considered. 

In Appendix A, we describe in detail the derivation of 
an effective master equation dp/dt — —i[H^p(t)] + ^ e ff(p)> 
where corresponds to the dipole force (10), and 



SW(p) = £ (Lfp(Lf)^ - {{LfVLfp - lp(Lf )tLf ) , 

n 

(36) 

accounts for the two possible decoherence channels. These 
are the so-called Raman and Rayleigh scattering of photons, 
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and are contained in the following effective jump operators 

+ yl (Q lt e i ( k >- r - ( »i ; ) + Q 2te i(k2T-<» 2 r)A | ;)(t | 5 

A V J (3?) 

Lf = ^ (^ 1 , t e i ( k '"- <Bir ) + Q 2it e i ( k2 "- <B2r )) |f>(t|+ 

+^ (n^C"*-^ +n 2 , ; e i ^--^)) |t>ai, 

where we have assumed that the scattering rate is much 
smaller than the laser detuning r«A. Note that the first 
term of each jump operator does not change the spin state (i.e. 
Rayleigh scattering), being responsible for pure dephasing. 
Conversely, the second term flips the spin state (i.e. Raman 
scattering), leading thus to damping of the spin populations. 

A conservative estimate of the photon scattering is to con- 
sider the individual effective rates (37), which scale as T e ff = 
r( | /A) 2 ~ |^l| (r/A). By using the parameters intro- 
duced in Sec. Ill, we find that r eff // eff « 10 3 (r/A). Hence, by 
considering a sufficiently large detuning, it is possible to keep 
the scattering rates below Tqr/Jqr < 10 -1 , so that they do not 
compromise the accuracy of the effective description (21). For 
instance, for typical decay rates T/2n ^1-10 MHz, one must 
consider detunings in the range A/27T ^10-100 GHz. 

(v) Spatial dependence of the laser-beam profile: The ef- 
fective spin couplings in Eq. (19) assume a Rabi frequency 
that is constant along the ion crystal. For large crystals, how- 
ever, the characteristic gaussian profile for the laser intensities 
may lead to weaker Rabi frequencies on the boundaries of the 
crystal. This effect will give rise to inhomogeneous spin-spin 
couplings that must be added to the inhomogeneities caused 
by the varying inter-ion distance in Coulomb crystals. Rather 
than considering these terms as an error, they can be seen as 
a gadget that makes the many-body problem even more inter- 
esting. In particular, they will lead to inhomogeneous critical 
points which may be responsible of interesting effects [46] . 

D. Efficient detection methods 

A crucial part of a QS is the ability to perform mea- 
surements that yield information about the Hamiltonian un- 
der study. For the frustrated quantum Ising ladders (22), 
a QS would start by preparing the so-called paramagnetic 
state |P) = <8>i|-»)i with |-^) f = (|t)/+ |^))/V2, which is the 
ground-state of the model in the absence of spin interactions. 
Such a separable state can be accurately prepared by means 
of optical pumping, followed by | -pulses globally addressed 
to the whole ion crystal [6] . This step is followed by the adi- 
abatic modification of the Hamiltonian parameters / e ff>^eff>^> 
in order to connect the paramagnet to other phases of the quan- 
tum Ising ladder (see Sees. V and VI below). 

A direct approach to the measurement of these phases is the 
so-called quantum state tomography, more precisely, the full 
determination of the state of the system [67]. However, this 
approach becomes highly inefficient for many-body systems 



due to the exponential growth of the composite Hilbert space, 
and alternative schemes must be studied. An interesting pos- 
sibility for state estimation are the methods based on matrix- 
product representations of the states [68]. Another alterna- 
tive that does not require full quantum state tomography is the 
measurement of order parameters characterizing the phases. 

One of the advantages of trapped-ion experiments with re- 
spect to other platforms is their ability to perform highly- 
accurate measurements at the single-particle level [6]. The 
technique of state-dependent fluorescence allows for the mea- 
surement of single and joint probability distributions of the 
electronic states pJ,pJJ. From the spatially resolved fluores- 
cence, it is possible to infer local expectation values, such as 
the magnetization (erf) = \{pj — 1), or two-body correlators 

(ofo-j) = \ [l -2(p} + p])+4P$]. As discussed in Sec. VI, 
these observables usually contain all the relevant information 
about the different phases. In particular, one could study the 
dependence of the correlator with the distance (48), or infer 
the magnetic structure factor S zz (q) = Y,ij{ G i G )) Qlq ^ l ~^ • 

Let us now comment on the possibility of recovering some 
of these magnitudes from global properties of the fluorescence 
spectrum. By measuring the probability to find a fraction of 
ft-ions in the excited state P^(n) [10], one obtains the total 
magnetization without single-site resolution m z = ^^-(of). 
To measure correlators, note that the fluorescence of an en- 
semble of emitters may carry information about their correla- 
tions [69]. In our case, the resonance fluorescence associated 
to the cycling transition depends on the collective properties 
of the ion ensemble [70] . In fact, the power spectrum in a par- 
ticular detection direction, f, is related to the structure factor 

S%(m) - £ e i ¥ f -( r ?- r °)((i + of )(i + o|)) oc s zz (q), (38) 

ij 

where X is the wavelength of the emitted light. Even if the ion 
spacing is much larger than the optical wavelength l z ^> A, one 
may compensate it by setting the photodetector almost orthog- 
onal to the plane defined by the ladder, and by using photode- 
tectors with a very good angular resolution. We finally note 
that such magnetic structure factors yield a lower bound on 
the entanglement without the need of state tomography [71]. 



V. SCOPE OF THE QUANTUM SIMULATOR 

Once the validity of the spin-ladder Hamiltonian (22) has 
been addressed by both analytic and numerical methods 
(Sees. II and III), and its experimental viability discussed 
(Sec. IV), we can now focus on the many -body models to 
be explored with the QS. We place a special emphasis on the 
range of collective phenomena that are not fully understood, 
or have not been addressed so far to the best of our knowledge. 
These would directly benefit from the advent of such a QS. 
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A. Ji-J 2 quantum Ising model 

The first many-body model that may be targeted with the 
proposed QS is the so-called axial next-to-nearest neigh- 
bor Ising model [78], supplemented by quantum fluctuations 
(see [34] and references therein). It has the Hamiltonian 

HM= J ^M+i+J2 1 £<y^ +2 -h , £af, (39) 

i i i 

which consists of nearest neighbor (J\ < 0) and next-to- 
nearest neighbor {J 2 > 0) couplings, and the transverse field h. 
In spite of the mild looking appearance of this Hamiltonian, it 
has a rich phase diagram with some features that are still con- 
troversial. We refer to this model as the J1-J2 quantum Ising 
model (/1-/2QIM), a paradigm of frustrated magnetism. 

The Hamiltonian of our QS (22) corresponds to the J\- 
/2QIM for the simplest possible geometry, namely, the two- 
leg zigzag ladder [18]. The indices of this ladder, 7= 1,2, 
4 = 1, • • • ,N/2 (Fig. 4), are mapped onto a one-dimensional 
chain by the relation i = 2(/ s — 1) + y (Fig. 5(a)). Then, the 
analogy with the /1-/2QIM follows directly provided that 

j = f\2 = Jeff COS 1 Jeff „ m 

1 J ™ |r?-f? +1 | 3 ' l^-^+al 3 

Note that due to the particular laser-beam arrangement pre- 
sented in Sec. II, it is possible to tailor the ratios J2/J1 and 
h/J\ experimentally. We emphasize that, even if a solid-state 
material is found to be described by such a model, a similar 
microscopic control of the couplings seems rather difficult to 
achieve. Therefore, the trapped-ion platform is ideal to ex- 
plore the different regions of the phase diagram [79]. Of par- 
ticular relevance is the region around the frustration point f c = 
Ji/\Ji\ = 1/^' which is characterized by a macroscopically- 
degenerate ground-state. Therefore, even a small amount of 
quantum fluctuations due to the transverse field may lift the 
classical degeneracy leading to a variety of magnetic phases. 

Such a rich phase diagram is studied numerically in Sec. VI, 
where we identify some additional features caused by the 
dipolar range of interactions present in the trapped-ion QS. 
Let us briefly mention that these long-range interactions intro- 
duce incompatible sources of frustration capable of stabilizing 
a new order that complements the ferromagnetic, dimerized 
antiferromagnetic, and floating phases that are also present 
in the short-range model (39). Hence, our QS will be of the 
utmost interest to explore the interplay between frustration, 
quantum fluctuations, and long-range interactions. Besides, 
the proposed QS shall be able to address some open ques- 
tions about the phase diagram that are still a subject of con- 
troversy [34], such as the extent of the floating phase and the 
existence of a multi-critical Lifshitz point. 



B. Dimensional crossover and quantum dimer models 

An ambitious enterprise is the understanding of the dimen- 
sional crossover from the two-leg quantum Ising ladder onto 



the two-dimensional (2D) triangular quantum Ising model 

#TQIM = Y< J2G m,n G m+l,n+ J l a m,n a m,n+l + J l 
mn 

— ^^hc m n , 

mn 

(41) 

where the spins are labelled according to the 2D Bravais lat- 
tice vectors r WjW = mai + n2L2 (see Fig. 4), and we consider 
anisotropic couplings J\ < and J2 > 0. The correspondence 
with our trapped-ion QS (22) is straightforward if one consid- 
ers that m labels the spins within each leg of the ladder cou- 
pled by J2 J, and n the different legs coupled by J\ o /. 
Note that, following recent experimental efforts [12, 80], the 
triangular QIM may also be realized with ions in Penning 
traps [81]. However, the study of the ladders and the crossover 
phenomena seems to be better suited to ions in Paul traps. 

The triangular classical Ising model can be considered as 
the backbone of frustrated magnetism [32]. Already in the ab- 
sence of quantum fluctuations, there are suggestive questions 
that deserve a careful consideration. For instance, the frustra- 
tion point of the two-leg zigzag ladder f c = 1/2 flows to the 
isotropic point f c = 1 in the 2D model. Moreover, the macro- 
scopic ground-state degeneracies of these models yield differ- 
ent ground- state entropies. We believe that it would be fasci- 
nating to explore these topics with trapped ions, which allow 
for the consecutive increase of the number legs (Figs. 1(c), (f)). 

The dimensional crossover is even more exotic when quan- 
tum fluctuations are included. Let us note that the dimer- 
ized antiferromagnet of the two-leg zigzag ladder corresponds 
to two possible classical dimer coverings of the ladder (see 
Sec. VI), where each dimer corresponds to a nearest-neighbor 
bond that is not satisfied due to the frustration. Since the 
ground- state tries to minimize the number of dimers, each site 
of the dual lattice belongs to only one dimer, and the covering 
consists of the dimer arrangement along the rungs of the lad- 
der. Note that each spin is connected to three satisfied bonds, 
and only one broken bond. The situation gets more interest- 
ing for the 2D quantum Ising model, since there, a single spin 
may be connected to the same number of satisfied and broken 
bonds. In this case, the transverse field can flip the spin and 
produce a resonating effect for neighboring dimers [37], pro- 
viding a beautiful connection to the so-called quantum dimer 
models (see [82] and references therein). 

Quantum dimer models were introduced [83] in the context 
oh high-temperature superconductivity. Here, these models 
provided a neat playground where to study the quantum spin 
liquid phases of the undoped cuprates, which were conjec- 
tured to play a key role in the onset of superconductivity upon 
doping [24]. However, they have evolved into an indepen- 
dent subject displaying exotic effects, such as topologically 
ordered phases and fractional excitations. Although originally 
introduced for Heisenberg antiferromagnets, where the dimers 
correspond to spin singlets, there is also a link to frustrated 
quantum Ising models by mapping the dimers to the broken 
magnetic bonds due to the frustration [37]. 

Their connection to the isotropic triangular quantum Ising 
model [37] has allowed to predict a quantum version of the 
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phenomenon of 'order by disorder' [84], which yields an or- 
dered phase out of the classically disordered frustration point 
f c = 1 when quantum fluctuations are switched on. In con- 
trast, by switching the transverse field on the frustration point 
f c = 1/2 of the two-leg zigzag ladder, one only obtains a 
disorder paramagnetic phase (see Sec. VI). Therefore, the 
trapped-ion QS offers a unique playground to study how the 
phenomenon of 'order by disorder' sets in as the number of 
legs is increased, and how the anisotropy and the long range 
of the interactions affects it. 



C. Quantum spin liquid phases 

In the models considered so far, strong quantum fluctua- 
tions trigger a phase transition connecting an ordered phase 
to the uninteresting disordered paramagnet. However, this be- 
havior does not exhaust all possibilities. Quantum fluctua- 
tions may be responsible of stabilizing more exotic phases that 
do not break any symmetry of the Hamiltonian, the so-called 
quantum spin liquids [85]. 

We have discussed how the two-leg zigzag ladder con- 
sisting of bond-sharing triangles allows for the QS of the 
paradigm of FQIM. However, this model only accounts for 
a disordered paramagnet. There exists another simple two-leg 
ladder which, although not so widely known, has been argued 
to provide the simplest instance of a frustration-induced quan- 
tum spin liquid [86] (see also [37]). This is the so-called saw- 
tooth quantum Ising model, which consists of corner- sharing 
triangles described by the following Hamiltonian 

#AA = J ill°i°Ui+ J 2L a 2i-i a 2i+i- h L°i- ( 42 > 

i i i 

Here, J\ < represents the coupling along the diagonal rungs, 
and fc > stands for the interactions along the lower leg of 
the ladder. In the isotropic frustration point f c = 1 , the ground 
state is found by minimizing the number of broken bonds 
in each triangle independently, which leads to a macroscopic 
ground-state degeneracy. The numerical results in [86], which 
are based on the exact diagonalization of small ladders and 
perturbative expansions, support the absence of any symme- 
try breaking as the transverse field is increased. This effect 
has been coined as 'disorder by disorder' , and would provide 
a testbed for a disordered quantum spin liquid state [37]. 

In order to perform a quantum simulation of this model, the 
geometry of a three-leg triangular ladder must be modified 
according to the method presented in Sec. IV (see Fig. 1(d)). 
Then, the model Hamiltonian (42) would follow directly from 
the QS Hamiltonian (22) with the usual identifications fc J, 
and J\ o /. In addition to the possibility of reaching larger 
spin ladders to test the conclusions of [86], the trapped-ion 
QS may explore of the effects of fc — fc/\J\ | 7^ fc and long- 
range interactions, hopefully leading to a rich phase diagram 
which, to the best of our knowledge, still remains unexplored. 

At this point, we should remark that quantum spin liq- 
uid phases can also be realized in (quasi) one-dimensional 
Heisenberg magnets [22, 26] . A much harder task is to find 
their higher-dimensional counterparts [85]. By using the same 



technique to modify the geometry of the trapped-ion ladder, 
one may construct three-leg ladders such as the one displayed 
in Fig. 1(e). Note that this amounts to a single stripe of the 
well-known two-dimensional (2D) Kagome lattice. Hence, 
our QS allows for the exploration of the dimensional crossover 
towards the Kagome quantum Ising model. In this case, the in- 
terplay between frustration and quantum fluctuations has been 
argued to give rise to a 2D quantum spin liquid phase [37], 
which could be targeted with the proposed QS. Besides, the 
study of the dimensional crossover and the long-range inter- 
actions is likely to introduce a variety of interesting effects. 

Before closing this section, let us remark that in addition to 
the aforementioned static phenomena, the trapped-ion QS is 
also capable of addressing dynamical many -body effects. In 
particular, the microscopic parameters of the above Hamilto- 
nians can be tuned dynamically across the different quantum 
phase transitions. The breaking of the adiabatic approxima- 
tion associated to a quantum critical point has been the subject 
of recent interest for different quantum systems (see e.g. [87]). 



VI. A DETAILED CASE: THE J x -fc QUANTUM ISING 
MODEL 



In this Section, we focus on the many-body physics of 
the /1-/2QIM. For the sake of completeness, we first review 
the properties of the short-range model, and then discuss the 
changes introduced by the dipolar range of the interactions, as 
realized in the trapped-ion quantum simulator. 



A. Short-range J\ -fc quantum Ising model 

The essence of the QS is captured by the idealized next-to- 
nearest neighbor QIM (39), rewritten here for convenience 

#M = -1*1 (lofof+i -/2l°f°f + 2 + sI°N . (43) 



where the ratios fc = Ji/\J\\ and g = h/\J\ \ can be experi- 
mentally tailored. The original J\-fc Ising model [78] con- 
siders competing ferromagnetic and antiferromagnetic inter- 
actions (see Fig. 6(b)), namely sign (72) = — sign(/i) = +1. 
According to the Toulouse- Villain criterion (20), this leads to 
a frustrated quantum spin model = — 1 > whereby the ratio 
/ controls the frustration, and g the quantum fluctuations. 

The classical J\-J 2 Ising chain, obtained by setting g — 
in the above Hamiltonian, can be solved exactly (see [88] 
and references therein). Such a solution yields two possible 
phases. For fc < \, one lies in a ferromagnetic (F) phase 
whose ground-state manifold is two-fold degenerate 

|F)espan{| t YV t ),| ; V--A)}, ( 44 > 

where the spins have been arranged according to the trapped- 
ion zigzag layout. Conversely, for fc > h, the ground-state is 
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Figure 10. Schematic phase diagram of the /1-/2QIM: Ferro- 
magnetic (F), dimerized anti-ferromagnetic (dAF), paramagnetic (P), 
modulated paramagnetic (mP), and floating (FP) phases. These states 
of matter are separated by different types of quantum phase transi- 
tions (Second Order, Kosterlitz-Thouless, and Pokrovsky-Talapov), 
and by the disorder line (dashed line). 



a dimerized antiferromagnet (dAF) with four- fold degeneracy 



|dAF)espan{| t y... t \) : 



; t... 4- 

it it 



1 4. t 4. t/' 1 4. t i/j • 



(45) 



Note that the degeneracy of the ferromagnetic manifold is re- 
lated to the global spin-flipping Z2 symmetry of the Hamil- 
tonian U = ®iOf, such that Z2 = {!,£/} is the smallest 
cyclic Abelian group. On the other hand, the doubling of 
the dimerized-antiferromagnet degeneracy is accidental (i.e. 
not related to symmetries). Such an accidental degeneracy 
becomes more important at the point f c = \ , where the frus- 
tration leads to a macroscopically degenerate ground-state. In 
fact, the degeneracy has been shown to scale exponentially 
with the number of spins d c (p N , where (p = |(1 + \/5) is 
the golden ratio [89]. This yields the hallmark of frustrated 
magnets, namely, a non- vanishing ground-state entropy. 

We are interested on the impact that quantum fluctuations 
may have on these degeneracies. For g ^> 1 ^> fa, the ground 
state corresponds to a single paramagnetic (P) state with all 
spins pointing towards the direction of the transverse field 



IP) = I- 



(46) 



where |— >•) = + \\)) / \fl. The situation gets more interest- 
ing for intermediate fields, whereby additional exotic phases 
and a variety of quantum phase transitions occur. Since the 
model is no longer integrable, the analysis of the full phase 
diagram has been a big challenge, requiring the combination 
of a variety of techniques. For instance, the mapping to a 
classical 2D model [90] yields a direct link to commensurate- 
incommensurate thermal phase transitions [88, 91]. Hence, it 
is possible to use the methods developed in this area in or- 
der to understand the magnetic phases of the quantum model, 
being numerical Monte Carlo [92] and free-fermion approxi- 
mations [93] two representative examples. Together with the 
more recent application of bosonization techniques [94] and 



Figure 1 1 . Conjectured phase diagram of the dipolar J\ -/2QIM: 

The longer-range dipolar Ising couplings, typical of Coulomb- 
mediated interactions, lead to the splitting of the multiphase point, 
and the appearance of a different ordered phase in between. 



numerical renormalization group methods [95], these stud- 
ies yield the rich phase diagram represented in Fig. 10. Note 
that in addition to the aforementioned phases, there is a mod- 
ulated paramagnetic phase (mP), and a highly-debated in- 
commensurate floating phase (FP). Since all these different 
phases meet at f c = \ , this macroscopically-degenerate point 
is also known as the multi-phase point. There are several crit- 
ical lines emerging from the multi-phase point, which give 
rise to quantum phase transitions of second order, Kosterlitz- 
Thouless [96], or Pokrovski-Talapov [97] type. 

Notwithstanding these big efforts, we emphasize that there 
still exists some controversy about the floating phase. In par- 
ticular, the extent of the floating phase is still a question of 
debate. Whereas some results point towards a FP that pro- 
longs towards fa ^> 1 , other treatments predict a finite region 
that terminates in a multi-critical Lifshitz point (see [34] and 
references therein). This makes a QS of the utmost interest to 
settle down these discrepancies. Moreover, as we discuss in 
detail below, the introduction of long-range interactions leads 
to additional open questions that have not been previously ad- 
dressed to the best of our knowledge. From our numerical 
survey, we conjecture that the dipolar-range of the couplings 
leads to the splitting of the multi-phase point and the appear- 
ance of an intermediate phase (see Fig. 11). 



B. Dipolar-range J1-J2 quantum Ising model 



Trapped ions are an ideal platform to test the effects of long- 
range Ising interactions. For non-frustrated systems, even 
if the model belongs to the same universality class as the 
nearest-neighbor case, these long-range interactions may shift 
the critical point, and favor long-distance quantum correla- 
tions [15]. For the frustrated systems under study, the effect 
of long-range interactions is expected to be more significant. 

By considering the dipolar range of the trapped-ion ladder 
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Figure 12. Incompatible sources of partial frustration: (a) The 

partial frustration for fa < \ in the short-range model can be mini- 
mized by one of the four possible dAF ground- states, (b) The intro- 
duction of longer-range couplings yields additional sources of frus- 
tration, which are incompatible with the dAF ground- state, and thus 
responsible of the splitting of the multi-phase point. 

Hamiltonian (22), the /1-/2QIM (43) must be modified to 

V i *Eodd (4?) 
i 8 e even i ' 

where we have introduced the ratios f§ = J§/\J\\, such that 
J§ account for the dipolar decay of the interactions. In this 
notation, the frustration ratio of the nearest-neighbor chain is 
f2=J2/\J\\, an d we can study the additional frustration com- 
ing from the dipolar range fa, fa.... We have analyzed numer- 
ically the phase diagram of such a long-range frustrated spin 
model by means of an optimized Lanczos algorithm that al- 
lows us to reach efficiently ladders with L = 24 spins. All the 
results showed below have been performed including the fa 
and fa terms of the dipolar tail, in addition to the competing 
nearest-neighbor (J\) and next-to-nearest-neighbor (J2) cou- 
plings. We have checked that the effect of including longer 
ranged interactions (i.e., those corresponding to fa and fa cou- 
plings) does not affect qualitatively the results. 

In order to distinguish the phases, we focus on the two-body 
correlators, which show the following scaling when \i — j\ ^> 1 
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where < mo < 1, £p, £mP stand for the correlation lengths of 
the paramagnets, 77 > characterizes the algebraic decay of 
the correlations in the gapless floating phase, and the different 
modulation parameters have been listed in Table I. 



Table I. Magnetic modulation parameters 
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According to these expressions, the F and d-AF dis- 
play long-range magnetic order with different periodicities, 
whereas the P and mP phases are disordered. Finally, the FP 
phase has quasi-long range order with a modulation parameter 
that flows with the ratios fa,g (hence the adjective floating), 
and is generally incommensurate with the underlying lattice. 

An observable capable of capturing the periodic modula- 
tions of the long-range ordered phases, and thus the features of 
the phase diagram, is the so-called magnetic structure factor. 
It is defined as the Fourier transform of the spin correlations 

S=(9)=I<of^>e i * (,w) , (49) 

ij 

with q G [0, 2k), and should attain a maximum at the different 
values shown in Table I for each of the phases. In order to 
evaluate the magnetic structure factor numerically, we have 
considered periodic boundary conditions, which shall capture 
the bulk properties in the center of the trapped-ion ladders. 

According to the scheme depicted in fig. 11, the main fea- 
tures of phase diagram of the dipolar /1-/2QIM can be divided 
into three different regions: the region where the ferromag- 
netic order prevails (fa < 0.5), the region with the dimerized- 
AF order (fa > 0.6), and the new intermediate region (roughly, 
fa G (0.5, 0.6)) that appears due to the competition of different 
long-ranged frustration mechanisms. In Figs. 12(a)-(b), we 
show how the dAF phase is destabilized by these competing 
mechanisms, which accounts schematically for the splitting of 
the multi-phase point. We will show below how the numerical 
evaluation of the structure factor supports this division. We 
note that the structure factor has been normalized to unity. 

In Fig. 13(a), we plot the order parameter S zz (^ ) corre- 
sponding to the dAF phase. It can be observed that both phase 
transitions, namely the discontinuous transition (Pokrovski- 
Talapov type) separating the dAF phase from the FP, and the 
continuous transition (Kosterlitz-Thouless type) between the 
FP and the mP phases, reflect themselves clearly in the struc- 
ture factor as two consecutive jumps with increasing trans- 
verse field. The fact that only one transition is observed for 
lower values of fa hint at the possibility that the multicriti- 
cal point has been shifted from g = to some finite value. 
Whether this shift survives for larger lattices, or corresponds 
to a finite- size effect, is an open question that cannot be ad- 
dressed due to the limitations of the diagonalization routine. 

The dependence of the structure factor with the lattice mo- 
mentum q can be seen in Fig. 13(b). Here, we represent the 
evolution of this magnitude as we increase g for a fixed value 
of fa. We have chosen some representatives for each phase: 
for low fields g, we are well into the dAF phase, whose peri- 
odicity is given by a single q = | component. Increasing the 
transverse field g, the state undergoes a transition to the FP, 
showing a (fa-, g-dependent) incommensurate modulation. 
Upon further increase of the field, we reach the mP phase. 
Note that for moderate g, the modulation of this phase is ap- 
parent in the peaked form of its structure factor. Increasing 
further the magnetic field polarizes the spins in the transversal 
direction yielding a flat vanishing structure factor. 

The same analysis has been carried out in the ferromagnetic 
region. In Fig. 14(a), we have plotted the 5^(0) component 
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Figure 13. Dipolar /i-/ 2 QIM (dAF region, L=24): (a) Antifer- 
romagnetic order parameter S zz (%/2). (b) Structure factor S zz (q) for 
.12 = 0.69. 



characteristic of a ferromagnet. In this case, however, only 
one clear phase transition arises as we increase the transverse 
field g for a fixed value of f 2 . Also in accordance with this 
observation, the modulation of the ground-state in Fig. 14(b), 
shows a shift between a pure ferromagnetic order for low val- 
ues of g, to a modulated paramagnetic one with increasing 
field, and finally to the completely polarized paramagnet for 
g ^> 1. This suggests the possibility that the unmodulated 
paramagnetic phase observed in the /1-/2QIM is not stable 
upon the effect of further dipolar couplings, or that its extent 
in the phase diagram is small enough to prevent to be accu- 
rately captured with finite- sized lattices. 

In the /1-/2QIM, there exists a one-point boundary between 
the ferromagnetic and dimerized-antiferromagnetic phases 
precisely located at the muticritical point (g = 0,/2 = 0.5). 
Interestingly enough, in the dipolar J\-J 2 QIM this is no longer 
true and both phases are separated by a finite region. In 
Fig. 15(a), we have plotted the structure factor S zz (q) for the 
lattice size L = 16 along a fixed value of f 2 in between the 
F- and dAF-phases. In this graph, a new type of modulation 
shows up in the form of two diferentiated peaks with momenta 
q\ = I and #2 = However, the precision in the determi- 
nation of such peaks is limited by the size of our lattice. We 
have carried out additional numerics for L = 24, which allow 
us to bound these peaks between q\ £ [f , f ] and q 2 G 7t\. 
We have checked that the relative amplitude of these peaks 
does not depend on the ratio between the intra- and inter-leg 
couplings. In Fig. 15(b), we plot the amplitude of these two 
characteristic peaks with increasing g. The shaded area repre- 
sents the range where both peaks coexist, i.e, the extent of the 
new ordered state. It would be very interesting to explore the 
precise origin of these modulations via the trapped-ion QS. 

We remark that the modulations q\ and q 2 coexist for low 
transverse fields. As we increase g, only the modulation with 
lower lattice momentum survives as a differentiated peak. In- 
deed, the momentum of this surviving modulation is compat- 
ible with that found in the mP phase for intermediate fields 
(compare the light-blue curve in Fig. 15(a) with the red ones in 
Figs. 14(b) and 13(b)). Whether this is a crossover or a quan- 
tum phase transition between the new conjectured phase and 
the mP, goes beyond the limitations of our numerical tools. 



Figure 14. Dipolar /1-/2QIM (F region, L=24): (a) Ferromagnetic 
order parameter S zz (0). (b) Structure factor S zz (q) for f 2 = 0.45 




Figure 15. Dipolar Ji-J 2 QIM (intermediate region): (a) Structure 
factor computed in a lattice of size L = 16. (b) Amplitude of the two 
differentiated peaks in the intermediate region. The accuracy of its 
location is conditioned by the size of the lattice: (solid lines) L = 24, 
qi = I and q 2 = (dashed lines) L = 16, q\ = |, q 2 = ^f. 



All the above results have been computed in ladders 

% = f (-1)*, fi = 0, z? = (i-\N)a, (50) 

with lattice parameters d = a . Note that a ladder composed 
of equilateral triangles would correspond to d = y/3a. The 
extent of the new intermediate phase in the dipolar J1-J2QIM 
is indeed strongly dependent on the geometry of the triangular 
plaquettes. A straightforward way of measuring the extent of 
this phase is computing the distance between the F and dAF 
phases along the line g = 0. In Fig. 16, we show the order 
parameters S zz (0) ,5^(^/4) and S zz (7i/2) for different values 
of the anisotropy ratio d/a, which correspond to the F phase, 
the new intermediate phase, and the dAF phase, respectively. 
From these graphs (qualitatively similar results are obtained 
with L = 24), it is apparent that the splitting of the multiphase 
point is strongly enhanced in anisotropic lattices. 



VII. CONCLUSIONS AND OUTLOOK 

Low-dimensional quantum Ising magnets are often repre- 
sentative models for a variety of emergent cooperative phe- 
nomena. Unfortunately, in contrast to their Heisenberg coun- 
terparts, the identification of materials accurately described by 
these models has turned out to be a much more difficult task. 
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Figure 16. Splitting of the multiphase critical point: order param- 
eters corresponding to the F (5^(0)) intermediate 0S zz (7z;/4))and dAF 
(5 zz (7r/2)) phases computed in a L = 16 lattice. Each figure corre- 
sponds to a fixed anisotropy ratio d/a. The homogeneous equilateral 
zigzag lattice corresponds to d/a = y/3 (figure at the bottom). 



Moreover, whenever the models include anisotropic interac- 
tions and tunable ladder structures, the prospects of realizing 
them diminish even further. A promising alternative to over- 
come such difficulties are the so-called quantum simulators. 

In this work, we have shown that cold ion crystals are 
promising candidates for the quantum simulation of a variety 
of quantum spin-ladder compounds. This avenue of research 
will allow for the study of collective phenomena due to the 
combination of frustration, quantum fluctuations, long-range 
interactions, and dimensional crossover. In particular, these 
ion crystals can be used to explore paradigmatic, yet contro- 
versial, models of quantum Ising ladders whereby analytical 
and numerical techniques seem to disagree. Moreover, they 
may also allow for the study of previously unexplored features 
of the models, such as the effects of long-range interactions. 

First, we have shown how to control the geometry of differ- 
ent self-assembled trapped-ion ladders. Based on this possi- 
bility, we have presented a protocol to tailor the anisotropy of 
the magnetic interaction mediated by the transverse phonons, 
which couple pairs of distant ions indirectly. When these ions 
belong to different legs of the ladder, it is possible to tune both 



the sign and the magnitude of the spin- spin couplings by ma- 
nipulating the laser-beam arrangement. The validity of this 
technique is supported by a detailed discussion of the possible 
sources of error in current ion-trap experiments, and by nu- 
merics showing an excellent agreement with our predictions. 

This tool opens a vast amount of possibilities for trapped- 
ion-based QS of cooperative magnetic phenomena. For in- 
stance, we have presented a thorough description of the QS for 
the cornerstone of frustrated quantum Ising magnets, the J1-J2 
quantum Ising model. Moreover, we have also discussed how 
the QS has the potential of realizing quantum Ising ladders 
with connections to the exotic quantum dimer models intro- 
duced in the context of high-temperature superconductivity, 
and address the intriguing dimensional crossover phenomena. 
Finally, we have also pointed out how this QS may yield a 
route towards the long- sought quantum spin liquid phases. 
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Appendix A: Spin-dependent dipole forces and dissipation 

In this Appendix, we present a detailed derivation of the 
effective laser-ion interaction [50, 51] in Eq. (10). By taking 
into account the spontaneous decay from the excited level (see 
Fig. 3), we can discuss the regime where the spin-dependent 
dipole forces arise, and also analyze the sources of error due to 
photon scattering in the experiments. We consider the master 
equation for the Lambda scheme in Fig. 3, namely 



dp 

df 



= -i[H + V,p(t)] + @(p(t)). 



(Al) 



The coherent part of the evolution is given by 



h = £e m |m)H, v = £ £ ^a ltS \r)(s\c- im "+n. c ., 

m=r,t4 /=l,2j=t,| 

(A2) 

where we have introduced the energies of the internal states 
£ n £|,£^, the Rabi frequencies of the transitions Q/ jjS , and the 
laser frequencies ft)/. The dissipator describing the sponta- 
neous decay from the excited state is of the Lindblad form 

®(P)= I {L nP Ll-\LlL nP -\ P LlL n ), (A3) 

w=l,2 

with the following jump operators L\ = y/T\\)(r\, and L2 = 
vf |t) (r\, where T is the spontaneous decay rate from the ex- 
cited state back to the spin manifold (see Fig. 3). Let us now 
define the detunings for all possible transitions 



Si s = £ r -£ s -(Di. 



(A4) 



As announced in the main text, when these detunings are 
much larger than the Rabi frequencies and the decay rate, 



21 



namely 5/ ?J ^> it is possible to adiabatically eliminate 

the excited state from the dynamics, and obtain an effective 
master equation within the spin manifold. We use the formal- 
ism introduced in [99], and obtain the master equation 

^ = -i[/W(0]+ 

+ Z( L fp(Lfr - HLf) f Lfp - \ P {LfVLf) , 

n 

(A5) 

where the Hamiltonian is // e ff = H s +// r + //d, such that 

H s = (e ; + A£ ; )||)(|| + (e r + A£ t )|t)(t| (A6) 
includes the following ac- Stark shifts 



Ae, 



\&Ls\ 2 Sls 



l=l,2*°l,s^ L t 



(A7) 



and we have introduced the sum of the two decay rates T t = 
2r. Additionally, we get the following two-photon stimulated 
Raman transitions 



if r = ^i^ / ,a-e i ^-^>+H.c, 



(A8) 



/,/' 



where we have introduced a = = (cr + )^, and the ef- 
fective Rabi frequencies for the different Raman transitions 



"(25 u +ir t )(25 / , t -ir t ) ; 



(A9) 



which involve the absorption of a photon from the beam I' = 
1 , 2 and posterior photon emission into the beam / = 1 , 2. The 
last term of the effective Hamiltonian is the one responsible 
for the spin-dependent dipole forces 

H d = L ^^l^^le^-^ + Hx, (A10) 

s=U 

where we have introduced the following Rabi frequencies 



(25 1 , 5 +ir t )(25 23 ,-ir t )- 



(All) 



Now, depending on the particular laser frequencies &>/, it is 
possible to select whether the laser-ion interaction leads to a 
stimulated Raman transition (i.e. (O^ := (0\ — Q>2 ~ &\- — £; =: 
C0o), or to a spin-dependent dipole force (i.e. co^ <C coo). The 
effects of the spin-dependent dipole force are more transparent 
by rewriting Eq. (A 10) as follows 

H d = |ft L Ie- ifl ^ + ^ L ^e- ic °L' +H.c, (A12) 

where we have introduced o z = |t)(t| — anc * 

+ ft L = (A13) 



Once the coherent part of the effective master equation has 
been derived, one must obtain the effective jump operators. In 
our setup, they can be expressed as follows 
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(A14) 



Note that the first term of each of the effective jump opera- 
tors corresponds to the so-called Rayleigh photon scattering, 
which takes place without modifying the internal spin state. In 
this formulation, it becomes clear why the Rayleigh scattering 
will only introduce dephasing when the amplitudes (i.e. terms 
between brackets) of each jump operator are different, as ob- 
served in recent experiments [77]. The second term of each 
jump operator corresponds to the Raman scattering, whereby 
the spin state is changed after the emission of the photon. 

In order to study the accuracy of this effective descrip- 
tion (A5), we confront it with the exact numerical integra- 
tion of the original master equation (Al). In Fig. 17(a), we 
compare both predictions for &^/e r = 0.1,£;/£ r = 0.05, and 
setting the laser parameters such that 5i ^/e r = 0.5, C0l = C0o, 

£2 2 ,t/ £ r = ^u/e r = °- 05 > and ^i,t = ^2,4. = 0- This set of 
parameters leads to the regime of stimulated two-photon tran- 
sitions, so that we expect to find periodic Rabi oscillations in 
the populations with a frequency ^ e ff/ £ r = 2.5 • 10 -3 when 
the initial state is \\jfo) = As shown in the upper panel 
of Fig. 17(a), these oscillations get damped due to the sponta- 
neous decay T/e r = 0.05, where a clear agreement of the exact 
and effective dynamics can be observed. In the lower panel, 
we represent the time-evolution of the coherences, which are 
damped due to the photon scattering. 

More interesting to our purposes are the results displayed 
in Fig. 17(b), where we have kept the same parameters as 
above, but set G)l = 10 -3 G^. This guarantees the absence 
of two-photon Raman transitions. By switching on the ad- 
ditional laser beams, such that Q^/^r = — ^it/^r = 0.05, 
the non- vanishing differential Rabi frequency in Eq. (A 13) 
should lead to the o z dipole force (A 12). Accordingly, we 
expect to find damped Rabi oscillations in the coherences for 
I Vb) = (| t) + I whereas the populations should only 

show an exponential damping. This agrees with the results 
displayed in Fig. 17(b), which also support the accuracy of 
the effective master equation in this regime. 

To recover the action of the lasers on the vibrational degrees 
of freedom, one should substitute £2/ s — >> £2/^e lk/ ' r in all the 
expressions above. Besides, in the regime of interest 

< coo < e r - e t , and T t <C S ijS « 5 ljt =: A, (A15) 

and summing over all the trapped ions, one recovers the dipole 
force term in Eq. (10) with the effective Rabi frequency = 
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Figure 17. Effective dissipative dynamics: (a) Population P±(t) = 

Tr{||)(||pO)}, and coherences {6 X ) = Tr{||)(t|e- ia} ° r p())} +c.c. 
in the regime of two-photon stimulated Raman transitions, (b) Same 
observables in the regime of spin-dependent dipole forces. 
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of strength £2l that only couples to the vibrational excitations, 
together with the ac- Stark shifts contained in (A7), must be 
compensated experimentally. Finally, the expression for the 
effective jump operators is the following 
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(A16) 

Since we are considering that \Q,^ S | <C A, it becomes clear that 
the effective scattering rates scale as T e ff = r(|£2/^|/A) 2 . 



Appendix B: Analysis of the ion micromotion 

In this Appendix, we take into account the micromotion by 
considering the radio-frequency (r.f.) trapping fields rather 
than the effective pseudo-potential in Eq. (1) [100]. More pre- 
cisely, the trapping potential is 

V t ({r 7 -}) = - f V cos(i2rf0 £ ( 1 + 4 tf-y})) 

if , X <B1) 

where Vb, C/o are the a.c. and d.c. potentials of the trap, ro,zo 
are the distances from the trap center to the electrodes and 
end-caps respectively, and K g < 1 is a geometric factor. To- 
gether with the Coulomb interaction, this trapping potential 
leads to a set of coupled Mathieu equations 

d 2 r ia £lr fl , „ e 2 x ,r ia -r ja 

-^- + -^[aa + 2q a cos^r* - - £ ^—^ = 0, 

(B2) 

where we have defined a x = 4eK g Uo/mzl^l 2 f = a y = —a z /2, 
and q x = — 2eVo /mr^Q^ = — q yj q z = 0. In the majority of the 
experimental setups, these parameters fulfill a a ,q 2 a ^ 1> so 
that one finds the following solution 

r ia = 4 (1 + ^acos(ftrff)) +Ar /a (0 (B3) 

where Ar^f) stand for the secular vibrations of the ions and 
frequency shifts due to micromotion [104, 105] that shall lead 
to the phonon Hamiltonian (9), and rf are the equilibrium po- 
sitions of the ion crystal. To lowest order in a a ,q 2 a , these are 
found from the static part of Eq. (B2), after solving 

f^(a a + ^)4- 2 l|^|=0, (B4) 

which is equivalent to Eq. (2) after identifying the effec- 
tive trapping frequencies co a = ^£l r f(a a + \ < fa) 1 ^ 1 • From 
Eq. (B3), one identifies two sources of radial micromotion. 
The second term corresponds to the micromotion associated 
to the small-amplitude secular vibrations of the ions, whereas 
the first term stands for an additional micromotion connected 
to the equilibrium positions lying off the trap axis. The for- 
mer can be minimized by laser cooling, whereas the latter is a 
driven motion that cannot be cooled, and is intrinsically linked 
to the planar structure of the ladder geometries (i.e. for linear 
ion chains, this micromotion can be compensated by align- 
ing the ions along the trap axis). Let us remark that for the 
regimes of interest, q x w 0.1-0.2, l z w 1-10 /im, this micro- 
motion can largely exceed that created by the secular oscil- 
lations, and must be thus considered as a potential source of 
error in our QS. The discussion below focuses on this type of 
micromotion r® a (t) = r® a (l + ^q a cos(Q lY ft)), which is usu- 
ally referred to as the excess micromotion [72]. 

Micromotion heating: One possible consequence of the mi- 
cromotion is undesired heating, either due to the r.f. field of 
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the ion trap, or to the additional lasers used for cooling. In 
particular, we focus on the transverse phonons since they are 
responsible for the spin-spin interaction. 

One term in the effect of micromotion can be understood 
from the interplay between the excess micromotion and those 
terms in the vibrational Hamiltonian (3) that do not conserve 
the number of vibrational excitations. By working in the local 
vibrational basis used in Sec. Ill, these terms amount to 

m) = ?I% V2 ^f (O^VJ^ + H-c, (B5) 

Z ij 

where the couplings i^ y (t) = (t) / (e 2 / 7|) correspond to 
those of Eq. (4) after taking into account the micromotion 



to(0 = 4 (1 + ^aCOs(ftrff)) . (B6) 



In order to neglect AH (t), which is responsible of the r.f. heat- 
ing, we expand Eq. (B5) to leading order in q a <C 1, and find 
that all the relevant terms can be neglected under a RWA if the 
following condition is fulfilled 



m « 



\2C0y dzQrfl £2 rf 



CO- 



CO- 



(B7) 



Since the ladder compounds present K y <C 1, and £2 r f ^> co z , 
the validity of the RWA is easily fulfilled. Thus the leading 
r.f. heating mechanism would be due to non linearities. This 
mechanism was analysed numerically in [101, 102], where 
scaling laws for the r.f. heating rates were predicted. In these 
studies, the heating rates originate from non linearities and 
thus scale with the initial temperature of the crystal. For the 
temperatures that apply for our scheme, the r.f. heating would 
be much smaller than the anomalous heating and thus could 
be neglected in our analysis. 

A different possibility is that of laser heating. Note that 
the QS requires laser cooling of the transverse phonon modes, 
although not necessarily to the ground-state. As emphasized 
in [72, 106], depending on the ratio of the r.f. frequency to 
the decay rate of the cooling transition, either a broadening of 
the transition or the appearance of multiple micromotion side- 
bands may occur, which can lead to undesired heating even 
when the laser frequency is tuned below the atomic resonance. 
To overcome this effect, one must carefully tune the laser fre- 
quency according to the regimes described in [72]. 

Micromotion contribution to the spin-dependent forces: An 
important question to address is whether the micromotion 
modifies the spin-dependent dipole force, as derived in Ap- 
pendix A. In the interaction picture, the Hamiltonian describ- 
ing the laser coupling (A2) for the whole ion crystal becomes 

V(t)=J^^\n)( Si \^^ (B8) 

l,s i 1 

after including the micromotion. Here, we have introduced 
%u = q x k[ • r 9/2, which represents the ratio of the radial excess 
micromotion (B3) to the wavelength of the laser radiation. 

To carry on with the analysis, we need to specify a partic- 
ular laser-beam arrangement (see the inset of Fig. 2(b)). We 



parametrize the laser wavevectors as follows 
k/ = — (cos a/e x + sin a/Cy), 



(B9) 



where a/ determines their angle with respect to the x-axis, 
and A sp is the wavelength of the ft 2 Si/2-^ 2 ^3/2 transition (see 
Fig. 3). Let us recall that in order to control the anisotropy of 
the spin interactions (19), Jfj f oc cos(0//), the corresponding 
angles must span the range fyj G [0,27r]. This implies that the 
laser-beam arrangement must fulfill 0C\ = — — Aa, where 
Aa <C | CCi | . One particular choice is the following 



Aa, Aa 

2 2d 



Asp «l, (BIO) 



where d is the inter-leg distance (see Fig. 4). With this choice, 
we find §2/ = 0, and % u = ±%q x , both fulfilling \% u \ <C 1. This 
property will allow us to truncate the following series 



,i^-cos(£V) 



(Bll) 



where J m (x) are the Bessel functions of the first kind. By 
substituting the series (Bll) in the laser-ion coupling (B8), 
we obtain a sum over all possible micromotion sidebands 

V(t) « ££^,,,7 m (^)ki)^|e ik ' r 'e i(5 '--" ,arf)r +H.c. 

l,s,i m 

(B12) 

According to Eq. (A 15), we can set A « 5/ jS , such that for a 
sufficiently large detuning £l rf /2n ^0.1 GHz <C A/2;r w 10 
GHz, we find that the micromotion sidebands only introduce 
a resonance for m ~ 10. However, the contribution of such 
terms is negligible since <C 1, and J m (^u) °* (£>u) m - As 
shown in [103], the correction to the leading term should be 
taken into account by adjusting the Rabi frequency. Let us 
also note that the remaining non-resonant sidebands can also 
be neglected in a RWA, since |^/^| <C A. Hence, we conclude 
that the validity of the spin-dependent dipole force derived in 
Appendix A is not compromised by the micromotion. 
To quantify the accuracy of this argument, we introduce 



£ m = max. 



{|(^(0)mic-(^(0)eff| 



t G 



<>'£]}' 



(B13) 



where (o x (t)} = Tr{||)(t|e~ ia) °'pM} + c.c. represents the 
coherences. In this expression, the effective evolution un- 
der the dipole force (cr x (r)) e ff, as represented in Fig. 17(b), 
is compared to the time evolution including the micromotion 
sidebands (<3" x (^)) m i c . Accordingly, £ m sets an upper bound to 
the error of the dipole force (10). In Fig. 18, we represent this 
error bound as a function of the relative micromotion ampli- 
tude We use the same parameters as in Appendix A, and 
set the r.f. frequency to £l r f/e r — 5 • 10 -3 , which is consistent 
with the constraint £l r f <C A. We observe that the micromo- 
tion sidebands only contribute with a small error (4-5%) for 
the regimes of interest %u = | q x , = (shaded region). 

Once this has been shown, we must consider the effects 
of the micromotion on the derivation of the spin-phonon 



24 




0.03 



0.25 



Figure 18. Micromotion error: Scaling of the micromotion error 
bound £ m with the ratio of the excess micromotion to the laser wave- 
length The shaded area represents the parameters of interest. 



term (12). Following a similar procedure, we obtain the mi- 
cromotion contributions to the Lamb-Dicke expansion con- 
necting the dipole force (10) to the spin-phonon coupling (12) 



H d = *f £ y m (^)e lk L r ?af/ p ({a„})e^ rf -'» L > + H.c. ) 

i,m,p 

(B14) 

where we have introduced £u = qx^-L • r f/2, and 

MM) = j } (^L^Jtk (^-^+£4^)^ . 

(B15) 

Note that the above expression (B 14) includes all the possible 
resonances between the secular and micromotion sidebands 
with the laser beatnote (Ol = C0\ — d>2- The leading-order term 
for small Lamb-Dicke parameter 7] nl _ <C 1 occurs for the sec- 
ular resonance m = 0, and p = 1 . This term leads directly to 
the desired spin-phonon coupling (12) provided that 



(B16) 



Considering the different orders of magnitude in the prob- 
lem col/271 « COy/2K «10 MHz <C ^ rf /27T « 0.1 GHz, the 
leading micromotion resonance would occur for the term 

&LTlZ ± Jm(%Li)j(~ pQ * -^rf-^)' with m = 1 and p « 9. Since 
these terms scale as (rf n ±) p with rf n ± <C 1, they get exponen- 
tially suppressed and can be thus neglected. Let us note that 
the remaining off-resonant sidebands can also be neglected 
via a RWA in the regime of interest <^ col <C £2 r f . There- 
fore, we conclude that the micromotion does not modify the 
spin-phonon coupling (12). 

Micromotion contribution to unwanted transitions: Equa- 
tion (B8) implicitly assumes that the dynamics due to the 
Raman-beam configuration can be accurately described by 
only three levels \U)> \ r i)- While this is always true for 
Zeeman ions, some special care must be taken for hyperfine 
ones, where the laser beams may excite some of the remaining 
states of the ground- state manifold due to the extra resonances 
introduced by the micromotion. Hence, the laser-ion Hamil- 



tonian in Eq. (B8) must be supplemented V + AV by 
AV(0 = £ ^|r I -)(fl I -|e^" r 'e i & OTn ^e i V + H.c., (B17) 

which includes all the additional states \ai) of the ground-state 
manifold with energies e a , and we have introduced the detun- 
ings Si A = e r — e a — COi. These detunings can be controlled by 
the Zeeman shifts of the ground- state manifold {\at)} caused 
by an external magnetic field. Note that the aforementioned 
unwanted transitions that take the state out of the spin sub- 
space s G {t, 1} follow from two-photon processes evolving 
like ^ 1 ,^*y m (^ 1 0/ m K^0*e- i( ^-^-(---^ rf >. There- 
fore, the Zeeman splittings 8 a , s = £ a — e s must be tuned to 



|fti, fl ftLI <|5 fl , 5 -CHL±Orf|, 



(B18) 



such that these unwanted transitions become highly off- 
resonant and can be neglected in a RWA. 



Appendix C: Thermal fluctuations and phonon heating 

In this appendix, we discuss an alternative derivation of the 
effective spin models based on the Heisenberg equation of 
motion, which shall allow us to predict the effects of finite 
temperatures for the results presented in Sec. III. The start- 
ing point is the time-independent spin-phonon Hamiltonian in 
Eq. (15), rewritten here for convenience 

% = + L(^» a f 4 + H.c), (Cl) 



where we have introduced &{ n = i%e lkL ' r ?7]^^^. In this 
picture, the evolution of the operators is given by the follow- 
ing system of coupled differential equations 



dt 

da n (0 
dt 

dof (0 
dt 



n 



(C2) 



= 0. 



Thanks to the last conserved quantity of(t) = of (0), we can 
integrate this system of equations exactly. The evolution of 
the phonon operators corresponds to that of a forced quantum 
harmonic oscillator, namely 



a n (t) = a j0) Q -^' + ^am (e"^'- l) . (C3) 



By substituting on the remaining equation, we find a homo- 
geneous linear differential equation for <7- + (f), which can be 
integrated exactly. We are interested in deriving an estimate 
for the scaling of the error at finite temperatures 



|<of)r-<of)r=o| 
\(°?)t=o\ 



(C4) 
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Figure 19. Relative thermal error: Scaling of the relative thermal 
error £j with the mean number of phonons in the center-of-mass 
mode n y . The analytical estimate £^ n in Eq. (C6) (red solid line) is 
compared to the exact numerical results £j x (doted lines) for different 
truncations of the phonon Hilbert space n t = 6, 7, 8. 



where T is the temperature determining the phonon Gibbs 
state p th = z-^-PZn^dan^ suc h that Z = Tr{e-PE« n « fl « fl «} 
is the partition function, and j3 = (&b^) _1 is expressed in 
terms of the Boltzman constant &b- By considering a sepa- 
rable initial state p(0) = | VsXVsl ® Pth> where | \j/ s ) is a pure 
spin state, we find the following expression for the relative 
thermal error 



£r= 1-Tr^p th e ^ 



(C5) 



This expectation value can be evaluated exactly to yield the 
following result 



e T = 1 



(l-cos^f)/*^ 



(C6) 



where the effects of the zero-point motion have been included 
in (of )r=o> and we have introduced the mean phonon num- 
bers for each of the vibrational modes = (a^a m ). It is 
interesting to note that, as argued for the so-called quantum 
phase gates [51], the error can be minimized by considering 
evolution times that are multiples of the detuning of the clos- 
est vibrational mode 5^* , namely tf = Inn/ 8^* , where n eZ. 

In order to check the validity of our derivation, we have 
confronted the prediction (C6) to the exact time evolution of 
the spin-phonon model in Eq. (CI). We use the parameters of 
Sec. Ill, namely, the trap frequency C0y/2n — 20 MHz, and the 
laser parameters r\ y = 0. 1 , C0l = 1 . 1 C0y, = 0. 15 1 S y \ /r\ y , and 
Cjc • = 0- In Fig- 19, we represent the exact results for the 
thermal error £^ x (dotted lines) obtained by the numerical in- 
tegration of the Liouville equation dp/dt = —i[H p + H&,p] up 
to tf — 7r/8/ e ff for the Hamiltonian (CI). Note that we truncate 
the Hilbert space of each vibrational mode to n t = 6 (yellow 
circles), n t = l (green squares), and n t = 8 (blue diamonds). 
These results are compared to the analytical estimate ef 1 (red 
solid line) in Eq. (C6), showing a remarkable agreement for a 
sufficiently large truncation of the phonon Hilbert space. 



Figure 20. Relative heating error: Scaling of the relative heating 
error £r h with the inverse of the heating rate T^ 1 (logarithmic scale). 



We now derive a useful expression for the scaling of the 
error in terms of experimental parameters, in the limit of 
\^in\ <^ 8^. We perform a Taylor expansion of Eq. (C6) 
for ic y « 1, and take into account the orthonormality prop- 
erties of the normal-mode displacements . For low- and 
high-temperatures, we have n„ w fiCl^ and n„ w 
respectively, where j3 = l/k^T. In both regimes, we find the 
following upper bound for the error 



St < % = max^je^} = 1 — e 



(C7) 



where we have introduced n y = (ficQy), which corresponds to 
the mean phonon number for the center-of-mass mode. Fi- 
nally, in the low-temperature regime, one finds 



4|a L pi|J 



(C8) 



By controlling the experimental parameters, this error term 
should be minimized for the QS. Note that this linear scaling 
of the relative error with the mean number of phonons coin- 
cides with the results shown in Fig. 19. 

Another important question to address is the error caused 
by heating mechanisms. The heating of a particular vibra- 
tional mode may be induced by a variety of factors, such as 
the combination of stray electric fields and fluctuating trap pa- 
rameters, elastic collisions with a background gas, fluctuating 
patch fields in the trap electrodes, or non-linear static electric 
fields (see [6] for details). Therefore, we shall not focus on a 
particular microscopic model, but use instead a phenomeno- 
logical master equation 



dp 

dt 



= -i[#p+fld,p(0] + ^(p(0), 



(C9) 



where the coherent dynamics is given by Eq. (CI), and the 
heating dissipator corresponds to 

^h(p) = £r h {a\pa n - \a n a\p - \pa n a\) . (C10) 
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Here, we have considered that the heating rate is equal for 
all the vibrational modes. In addition to Eq. (C6), the heating 
mechanism provides another source of error. In order to single 
out such a contribution, we consider an initially ground-state- 
cooled ion crystal, and integrate numerically the above master 
equation. To estimate the relative error caused by the heating 
mechanism, we evaluate the following figure of merit 



l(of>r h -(of>r h =o| 



(Cll) 



l(^)r h =o| ' 
By considering a small timescale tf <C (Th) -1 , the evolution 



of the mean number of phonons due to the dissipator (CIO) 
yields a simple linear heating n„(t) = {a\a n ) « T\J. We set 
tf= \nj (8/ e ff), in such a way that n^(t) <C n t = 2. The re- 
maining parameters are the same as above. In Fig. 20, we 
represent the relative error (Cll) as a function of the inverse 
heating rate T^ 1 , which sets the timescale for the creation of 
one vibrational excitation. As can be seen in this figure, heat- 
ing times above 5 ms/ phonon, only have a small contribution 
(1%) to the overall error of the QS. 
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